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It  will  also  reduce  time  in  the  basin  and  more  accurately  represent  a  full-scale  seaway  where 
ship  may  be  at  risk. 

To  investigate  the  feasibility  of  producing  grouped  extreme  waves  in  an  experimental 
basin,  experiments  were  conducted  in  August  and  November  2007  in  the  Maneuvering  and  SeakeepiiB 
(MASK)  basin  at  the  Naval  Surface  Warfare  Center,  Carderock  Division  (NSWCCD) .  The  fir® 
experiment  (Phase  I)  consisted  of  producing  several  combinations  of  finite  regular  waves  with 
varying  parameters,  including  amplitude,  frequency,  and  signal  duration.  Superposition  of  the 
finite  regular  waves  produced  repeatable  wave  groups  at  a  specific  location  in  the  MASK.  T 
second  experiment  (Phase  II)  consisted  of  two  parts.  The  first  part  applied  the  wave-pack 
method  to  produce  single  large-amplitude  waves,  based  on  the  technique  developed  by  Clauss  and 
Bergmann  (1986)  and  Clauss  and  Kuehnlein  (1994,  1995).  The  second  part  employed  t 

achievements  from  Phase  I  to  embed  grouped  extreme  waves,  obtained  through  finite  regular  wa 
superposition  in  two  scaled  irregular  seaways:  a  Bretschneider  sea  state  8  spectrum  and  a 
Hurricane  Camille  spectrum. 

The  maximum  calibrated  attained  wave  height,  normalized  by  a  typical  ship  model  leng 
(h/L)  was  0.205,  compared  to  0.1  for  the  typical  largest  regular  wave  height  to  ship  mod 
length  ratio  used  for  testing,  or  double  the  maximum  wave  height  normally  achieved  in  the  MASK. 
The  maximum  wave  steepness  (h/A)  observed  was  approximately  1/7,  approaching  the  theoretic 
limit  for  a  non-breaking  wave,  compared  to  a  1/10  wave  steepness  which  is  typically  the  maxim 
for  regular  wave  testing.  The  largest  estimated  full-scale  wave  height  observed  was  37.5  m  (1 
ft.)  at  a  scale  ratio  of  46.6.  Producing  repeatable  wave  groups  with  the  desire 

characteristics  at  a  fixed  location  is  a  non-trivial  task,  however,  this  study  demonstrated  t 
feasibility  of  producing  groups  of  asymmetric  extreme  waves  in  the  MASK.  Future  work  is  al 
discussed . 
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Abstract 


Experiments  were  performed  to  develop  the  capability  to  generate  groups  of  large- 
amplitude  waves  in  irregular  seas  in  a  model  basin.  Secondary  objectives  included  demonstrating 
the  ability  to  generate  single  extreme  waves  and  also  the  application  of  experimental 
measurement  techniques  to  produce  point  and  wave-field  topology  data.  The  process  was 
intended  to  be  deterministic  in  nature,  such  that  a  large  wave  group  or  single  large  wa\e  will 
occur  in  the  model  basin  at  a  predictable  and  repeatable  location  and  time.  Proving  the  feasibility 
of  generating  asymmetric  large-amplitude  wave  groups  in  an  experimental  basin  is  the  first  step 
in  the  development  of  an  experimental  test  technique  that  ensures  a  model  will  be  exposed  to 
multiple  realistic  extreme  wave  events  during  a  test  run.  This  technique  will  remove  the 
“randomness"  associated  with  model  experiments  in  irregular  waves  from  the  proeess  to  e\aluate 
ship  performance  in  severe  sea  conditions.  It  will  also  reduce  time  in  the  basin  and  more 
accurately  represent  a  tull-scale  seaway  where  a  ship  may  be  at  risk. 

To  investigate  the  feasibility  of  producing  grouped  extreme  waves  in  an  experimental 
basin,  experiments  were  conducted  in  August  and  November  2007  in  the  Maneuvering  and 
Seakeeping  (MASK)  basin  at  the  Naval  Surface  Warfare  Center,  Cardcrock  Division 
(NSWCCD).  The  first  experiment  (Phase  I)  consisted  of  producing  several  combinations  of  finite 
regular  waves  with  varying  parameters,  including  amplitude,  frequency,  and  signal  duration. 
Superposition  of  these  finite  regular  waves  produced  repeatable  wave  groups  at  a  specific 
location  in  the  MASK.  Ihe  second  experiment  (Phase  II)  consisted  of  two  parts.  The  first  part 
applied  the  wave-packet  method  to  produce  single  large-amplitude  waves,  based  on  the 
technique  developed  by  Clauss  and  Bergmann  (1986)  and  Clauss  and  Kuehnlein  (1994,  1995). 
The  second  part  employed  the  achievements  from  Phase  I  to  embed  grouped  extreme  waves, 
obtained  through  finite  regular  wave  superposition  in  two  sealed  irregular  seaways:  a 
Bretsehneider  sea  state  8  spectrum  and  a  Hurricane  Camille  spectrum. 

The  maximum  calibrated  attained  wave  height,  normalized  by  a  typical  ship  model  length 
(h/L)  was  0.205,  compared  to  0. 1  for  the  typical  largest  regular  wave  height  to  ship  model  length 
ratio  used  for  testing,  or  double  the  maximum  wave  height  normally  achieved  in  the  MASK.  The 
maximum  wave  steepness  (h/X)  observed  was  approximately  1/7,  approaching  the  theoretical 
limit  for  a  non-breaking  wave,  compared  to  a  1/10  wave  steepness  which  is  typically  the 
maximum  for  regular  wave  testing.  The  largest  estimated  full-scale  wave  height  observed  was 
37.5  m  (123  ft.)  at  a  scale  ratio  of  46.6.  Producing  repeatable  wave  groups  with  the  desired 
characteristics  at  a  fixed  location  is  a  non-trivial  task,  however,  this  study  demonstrated  the 
feasibility  of  producing  groups  of  asymmetric  extreme  waves  in  the  MASK.  Future  work  is  also 
discussed. 
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Nomenclature 

h  wave  height 

H,  significant  wave  height 

hfs  estimated  full-scale  wave  height  for  notional  model 
/,  notional  scale  model  length 

X  wavelength 

T  wave  period 

k  wave  number 

s  wave  steepness 

cg  group  velocity 

f  frequency 

a  Rayleigh  distribution  scale  parameter 

xa  concentration  position 

ta  concentration  time 

sn  normalized  spectrum 

Hn  amplitude  of  wave-maker  piston  motion 
Cyc„  number  of  cycles  in  finite  regular  wave 
ua  Type  A  standard  uncertainty 

U  expanded  uncertainty 

kf  coverage  factor 

r  correlation  coefficient 


2 


Introduction  and  Background 


This  teehnieal  report  details  progress  made  in  FY07  and  FY08  for  a  sub-task  of  the 
Surfaee  Ship  Dynamics:  Motions  and  Maneuvering  projeet,  sponsored  by  the  Office  of  Naval 
Research.  The  sub-task.  Definition  of  Towing  Basin  Environment,  is  part  of  a  combined  effort  at 
NSWCCD  to  improve  predictions  and  measurement  of  ship  motions  in  waves  and  assess  the 
dynamic  stability  and  seakeeping  perfonnanee  of  naval  ships. 

The  advent  of  novel  hull  form  geometries  has  demonstrated  a  need  for  improved 
simulation  capabilities  and  experimental  techniques.  These  are  necessary  to  assess  stability 
perfonnanee,  to  increase  understanding  of  the  principle  physics  governing  dynamic  stability 
events,  and  to  provide  correlation  data  for  the  development  of  simulation  tools. 

Beeause  of  the  potential  risk  to  oeean-going  vessels  and  platforms  eaused  by  extreme 
seas,  an  increase  in  the  understanding  of  the  characteristics  and  occurrence  of  large  waves  is 
important.  Large  amplitude  waves,  either  singularly  or  in  groups,  can  present  a  serious  stability 
risk  to  a  ship,  leading  to  payload  and  hull  damage,  personnel  injury,  or  in  the  worst-ease 
scenario,  eapsize  and  loss  of  life  and  the  vessel. 

Large-amplitude  waves  have  been  the  subject  of  eye-witness  accounts  at  sea.  A  report  by 
Admiral  Robert  Fitz-Roy,  published  in  1839,  recounted  the  HMS  Thetis  traveling  between  two 
waves,  believed  to  be  greater  than  18.3  m  (60  ft)  wave  height,  beeause  the  wave  crests  were 
taller  than  her  masts  and  the  sails  were  slaek  despite  heavy  winds.  Other  large  waves  believed  to 
be  reliably  reported  inelude  a  34  m  (112  ft)  wave  reported  by  the  USS  Rampo  in  the  North 
Paeifie,  in  1933,  and  a  27.5-30.5  m  (90-100  ft)  wave  that  struek  an  oil  rig,  near  Vancouver 
island,  in  1968.  Additional  reports  have  been  made  of  ships  sailing  into  “holes  in  the  sea,"  made 
by  the  troughs  of  large  waves  (Draper.  1971;  Buckley  2005). 

These  large  waves  are  known  eommonly  by  several  names,  such  as  freak,  rogue,  or 
extreme  waves.  For  the  purpose  of  this  investigation,  they  will  be  referred  to  as  extreme  waves. 
They  typieally  are  non-linear  and  asymmetrie,  featuring  deep  troughs,  either  preceding  or 
following  steep  wave  erests. 

Research  has  indicated  the  relatively  common  occurrence  of  these  extreme  waves,  much 
more  often  than  was  predieted  by  conventional  Rayleigh  distribution  prediction  models 
(Rosenthal  and  Lehner  2008).  A  European  Union  collaboration  projeet.  Max  Wave  (Rosenthal 
and  Lehner,  2004,  2008),  began  to  investigate  extreme  waves  in  2000.  During  a  three-week 
period  in  2001,  European  Spaee  Ageney  (ESA)  satellites  recorded  data  from  an  imaged  area  of 
1 .5  million  square  kilometers.  Wave  heights  of  25  m,  or  larger,  were  observed  at  an  equivalent 
rate  of  one  in  150,000  square  kilometers  in  the  three  week  period. 

Approaches  for  estimation  of  the  probability  of  oeeurrenee  of  single  extreme  waves  in  an 
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experimental  basin  or  ocean  environment,  at  either  a  fixed  point  or  within  a  wave-field,  have  also 
been  investigated  (Lopatoukhin,  et  al.  2004;  Fedele,  2006a,  2006b).  Guedes  Soares  and  Pascoal 
(2005)  and  Petrova,  et  al.  (2007)  showed  the  inadequacy  of  design  wave  models  based  on  linear 
assumptions  and  the  difficulty  of  second-order  hydrodynamic  models  to  describe  accurately 
large  ocean  wave  characteristics,  as  determined  from  field  measurements. 

Extreme  waves  are  typically  characterized  by  wave  heights  of  at  least  two  times  the 
significant  wave  height.  Although  this  criteria  is  widely  used,  it  is  not  a  universally  accepted 
definition.  A  more  general  definition,  proposed  by  Johannessen  and  Swan  (2003)  states  that  an 
extreme  wave  can  be  characterized  as  a  wave  with  a  crest  elevation  significantly  larger  than 
considered  by  the  statistics  of  the  underlying  frequency  spectrum.  These  waves  of  large- 
amplitude,  either  single  or  grouped,  have  been  reported  in  all  the  oceans  of  the  world  (Rosenthal 
and  Lehner,  2008). 

Measurements  were  also  made  over  a  six-year  period,  from  1098-2003,  on  a  platform  in 
the  southern  Indian  Ocean  (Liu  and  MacHutchon,  2008).  This  region,  off  the  east  coast  of  South 
Africa,  is  a  major  shipping  route  and  is  known  to  produce  extreme  waves  because  of  the  mix  of 
the  Agulhas  Current  and  storm  waves  (Sverdrup,  et  al.,  2003).  These  measurements  have 
provided  a  preliminary  data  set  for  a  possible  class  of  extreme  waves,  characterized  by  wave 
heights  of  at  least  four  times  the  developed  significant  wave  height  (Liu  and  MacHutchon  2008). 

Extreme  waves  have  been  shown  to  occur  in  all  sea  states  (Bitner-Gregersen  and  Hagen 
2004),  although  they  are  of  greater  concern  to  ship  safety  in  higher  sea  states.  Some  research  has 
indicated  the  possibility  that  extreme  waves  may  be  more  likely  to  form  in  seas  dominated  by 
swell,  than  in  a  locally  wind  generated  sea  (Gibson  and  Swan,  2007).  Due  tc  the  complex 
interaction  of  swell  and  wind-generated  waves,  the  exact  environmental  mechanisms  for 
producing  extreme  waves  are  still  not  well  understood.  Possible  causes  of  extreme  waves 
forming  singularly  or  in  groups  include  well-developed  storm  systems,  intersection  of  wave 
systems  from  different  storm  centers,  and  storm  situations  where  the  storm  and  wave  group 
velocities  are  similar  (Rosenthal  and  Lehner  2008). 

The  research  to  date  has  focused  primarily  on  detection,  observation,  measurement,  and 
reproduction,  either  experimentally  or  numerically,  of  single  extreme  waves.  The  existence  of 
these  singular  large-amplitude  waves,  has  only  been  confirmed  by  measurements  within  the  last 
two  decades,  such  as  the  famous  “New  Year’s  Wave”  recorded  from  the  Draupner  platform  in 
the  North  Sea  (Figure  1).  In  2000.  in  the  Rockall  Trough,  west  of  Scotland,  a  British 
oceanographic  research  vessel  recorded  waves  with  a  significant  wave  height  of  18.5  m  (61  ft) 
and  individual  waves  up  to  29. 1  m  (95  ft).  These  waves  are  the  largest  ever  measured  in  the  open 
ocean  (Holliday,  et  al.,  2006). 
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Figure  1.  The  Draupner  Platform  wave,  January  1,  1995,  Time-history  of  the  wave 
amplitude,  in  meters  (Rosenthal  and  Lehner,  2004) 

During  Hurricane  Ivan,  measurements  made  by  Naval  Research  Laboratory  instruments 
on  the  floor  of  the  continental  shelf  in  the  Gulf  of  Mexico,  recorded  significant  wave  heights  of 
17.9  m  (58.8  ft)  and  individual  maximum  wave  heights  of  27.7  m  (91  ft).  NRL  researchers  ha\e 
stated  that,  although  not  recorded,  additional  analysis  indicated  significant  wave  heights  may 
have  exceeded  21  m  (69  ft),  with  individual  wave  heights  exceeding  40  m  (132  ft)  during  the 
storm  (Wang,  et  al.,  2005).  However,  within  hurricane  regions,  small  data  samples  have  been 
shown  to  underestimate  wave  heights  for  extreme  waves  (Jonathan  and  Ewans,  2007),  resulting 
in  the  possibility  that  current  measurements  may  not  even  provide  an  accurate  indication  of  the 
upper  limit  of  extreme  wave  heights. 

Additional  accounts  of  ships  experiencing  groups  of  exteme  waves,  such  as  the  “Three 
Sisters,”  have  been  reported  (Buckley,  2005;  Smith,  2006).  A  wave  group  is  defined  as  a  series 
of  waves,  with  wave  heights  larger  than  a  specified  threshold,  and  with  approximately  equal 
periods  (Masson  and  Chandler,  1993;  Ochi,  1998).  A  2005  U.  S.  National  Transportation  Safety 
Board  investigation  conducted  on  a  maritime  accident  in  the  Atlantic  Ocean,  off  the  coast  of 
Georgia  involving  the  Norwegian  Dawn — indicated  that  the  vessel  had  encountered  three  large 
waves  in  succession  (NTSB,  2005).  Experimental  generation  of  groups  of  extreme  waves,  has 
only  focused  on  a  symmetric  three  wave  group  of  heights,  Hs,  2HS,  Hs  (Clauss,  2002a).  While 
useful  for  determining  critical  motions,  the  symmetric  wave  group  does  not  represent  many  of 
the  characteristics  of  extreme  wave  groups  observed  in  the  ocean  environment. 

To  model  extreme  ship  motions  experimentally,  a  large  test  matrix  of  runs  must  be 
performed.  Model  experiments  conducted  in  random  waves,  to  assess  severe  ship  motions, 
dynamic  stability,  slamming,  and  ultimate  strength;  require  long  test  am  times  to  ensure  that  the 
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low  probability  extreme  wave  events,  which  still  occur  in  the  seaway,  have  been  encountered. 

The  first  approach  to  reduce  the  number  of  required  tests  for  regular  wave  seakeeping 
was  the  transient  wave  technique,  developed  analytically  by  Davis  and  Zamick  (1964).  At  the 
David  Taylor  Model  Basin,  Davis  and  Zamick,  and  Gersten  and  Johnson  (1969)  applied  the 
transient  wave  technique  to  regular  wave  model  experiments  for  heave  and  pitch,  at  zero  and 
forward  speed.  These  tests  demonstrated  a  potential  reduction  by  an  order  of  magnitude  of  the 
total  necessary  testing  time.  The  transient  wave  technique  was  also  applied  to  model  testing  in 
Japan  in  the  mid-1970s  (Takezawa  and  Takekawa,  1976;  Takezawa  and  Hirayama,  1976). 

Clauss  and  Bergmann  (1986),  Clauss  and  Kuehnlein  (1994,  1995),  and  Matos,  et  al. 
(2005)  used  a  transient  wave  technique  with  Gaussian  wave-packets  to  excite  model  ships  and 
offshore  structures  in  an  experimental  basin.  Further  revisions  to  this  technique  included  using 
nonlinear  transient  wave  trains  and  modified  wave  celerity  to  generate  extreme  waves  (Clauss, 
1999).  Transient  waves  have  even  been  proposed  for  use  as  part  of  a  wave-maker  calibration 
procedure  (Masterton  and  Swan,  2008).  Model  tests  have  been  performed  with  large  transient 
waves,  embedded  in  both  regular  and  random  wave  trains.  These  waves  were  calculated  both 
linearly  and  nonlinearly,  with  empirically-based  terms  for  the  particle  orbital  motion  and 
shallow- water  effects  (Clauss  and  Hennig,  2002).  Experimentally,  single  extreme  waves  have 
been  shown  to  have  larger  asymmetry  than  full-scale  storm  waves,  but  with  a  profile  more 
similar  to  real  extreme  waves  than  second-order  numerical  simulations  (Antao  and  Guedes 
Soares,  2008). 

To  generate  desired  deterministic  wave  sequences  in  the  experimental  basin,  techniques 
where  singular  extreme  waves  were  embedded  in  irregular  seas  with  a  linear  wave  theory  “first 
approach,”  and  then  optimized  using  a  fully  nonlinear  approach  have  been  investigated  (Clauss, 
2000;  Clauss,  2002a;  Clauss  and  Schmittner,  2005).  Experiments  using  a  modified  “New  Year 
Wave”  have  been  conducted  to  assess  both  motions  and  structural  response  for  floating  offshore 
structures  (Clauss,  et  al.,  2008).  Model  tests,  to  induce  extreme  roll  and  capsize  for  a  ship,  must 
consider  wave  characteristics  such  as  wave  height  and  steepness,  wave  groupiness,  and  the 
velocity  and  direction  of  wave  propagation  (Clauss  and  Hennig,  2004). 

Envelope  theory  (Longuet-Higgins,  1957),  used  to  mathematically  describe  the  ocean 
wave  environment,  assumes  the  wave  spectra  is  Gaussian  and  narrow-banded.  Kimura  (1980) 
proposed  describing  a  series  of  waves  above  a  specified  threshold,  using  a  Markov  chain  to 
obtain  the  probability  distribution  of  wave  groups  in  a  random  sea.  Additional  revisions  were 
proposed  by  Longuet-Higgins  (1984)  to  simplify  the  application  of  Markov  theory  to  describe 
the  statistics  of  steep  wave  groups,  and  by  Battjes  and  Van  Vledder  (1984)  to  link  the  wave 
group  properties  to  the  energy  spectrum.  Neither  spectral  method  is  able  to  fully  describe  the 
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groupincss  characteristics  of  a  given  wave-field  (Masson  and  Chandler,  1993).  However,  recent 
analysis  of  environmental  data  from  the  North  Sea  has  shown  the  ability  of  a  Markov  chain 
model  to  better  predict  wave  group  properties  than  a  Rayleigh  distribution  model,  except  for 
large  wave  heights  (Stansell,  et  al.,  2002).  Because  of  the  Gaussian  assumption,  application  of 
both  envelope  and  Markov  theory  may  be  limited  for  realistic  extreme  waves  and  extreme  wave 
groups.  Directionality  has  also  been  shown  to  be  an  important  factor  for  accurately  modeling 
extreme  waves  observed  in  the  ocean  (Johanncssen  and  Swan,  2001;  Gibson  and  Swan,  2007). 
Other  analytical  methods,  based  on  the  nonlinear  Schroedinger  equation  (Osborne,  2001,  2006; 
Slunyacv,  et  al.,  2005;  Grue,  2002)  and  the  Kadomtsev-Petviashvili  equation  (Porubov,  et  al., 
2005),  have  been  used  both  to  model  single  large-amplitude  waves  in  an  experimental  basin  and 
to  numerically  recreate  full-scale  ocean  wave  measurements. 

“Numerical  wave  tanks”  have  been  developed  for  the  deterministic  analysis  of  ocean 
structure  behavior  in  simulated  single  extreme  waves  (Clauss,  et  al.,  2005;  Ning,  et  al.,  2008). 
One  numerical  tank  applies  potential  flow  codes  for  fast  computations  and  RANS  codes  to  model 
breaking  waves  and  fluid-structure  interaction  (Clauss,  et  al.,  2005).  Some  agreement  exists 
between  computations  and  experiments,  with  potential  flow  and  RANS  used  either  separately,  or 
coupled.  However,  wave  heights  are  shown  to  be  generally  over-predicted  in  a  numerical  wave 
tank.  Extreme  waves  simulated  in  numerical  tanks  have  also  been  applied  to  predict  structural 
loading  on  ships  (Guedes  Soares,  et  al.,  2008)  and  offshore  structures  (Clauss,  2002b).  However, 
loads  predictions  have  been  shown  to  vary  widely,  depending  on  the  wave  kinematics  model 
employed  for  the  evaluation  (Stansberg,  et  al.,  2008),  models  such  as  second-order  (Bitner- 
Grcgcrsen  and  Hagen,  2004),  Grue’s  method  (Grue,  ct  al.,  2003),  or  the  well-known  Wheeler 
stretching  method  (Wheeler,  1969).  Some  comparisons  between  numerical  predictions  and 
extreme  wave  model  experiments  were  performed  by  Hennig,  et  al.  (2006). 

Tools  for  ship  design  have  been  developed  where  waves  or  wave  groups  arc  used  to 
induce  a  specific  ship  motion  response.  This  approach  was  first  applied  by  Tikka  and  Paullmg 
(1990),  using  wave  groups  to  induce  large  roll  excitations  caused  by  parametric  roll.  Additional 
studies  of  applications  of  wave  groups  to  parametric  roll  response  have  been  made  by 
Boukhanovsky  and  Degtyarev  (1996).  Alford  has  used  a  design  wave  train  method  to  produce  a 
desired  motion  response  (Alford,  et  al.,  2006;  Alford,  2008).  The  method  has  also  shown  that 
wave  grouping  and  directionality  are  significant  factors  in  the  realization  of  worst-case  ship 
motions  (Alford,  et  al.,  2007).  To  determine  particular  severe  instabilities  for  a  ship,  Thcmclis 
and  Spyrou  (2007)  applied  critical  wave  groups  deterministically  by  predicting  the  required 
critical  wave  groups  to  induce  the  instability,  for  a  given  ship.  Then  the  probability  of 
encountering  one  of  these  critical  wave  groups,  for  a  given  route  and  time,  was  computed.  Using 
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this  critical  wave  group  method,  instabilities  were  assessed  including  beam-seas  resonance, 
parametric  roll,  and  pure  loss  of  stability. 

All  wave  conditions  do  not  have  the  same  probability  of  occurrence.  To  predict  rare 
response  events,  wave  groups  can  present  a  scenario  of  higher  probability  for  extreme  response 
than  a  single  large-amplitude  wave.  For  both  experiments  and  simulations,  realistic  groups  of 
large-amplitude  waves  can  be  applied  to  overcome  the  “problem  of  rarity”  (Belenky,  et  al.,  2008) 
by  inducing  realistic  severe  conditions  for  large  roll  motions  or  stability  failure,  at  a  known  time 
and  location. 

During  large-scale  model  testing  at  Patuxent  River  Naval  Air  Station  (NAS)  in  the  fall  of 
2006  (Carrico,  2007),  wave  time-histories  taken  from  buoys  deployed  in  the  Chesapeake  Bay 
testing  site  showed  waves  with  characteristic  groupiness  (Figures  2-5).  Both  singular  and 
grouped  extreme  waves  were  present  in  some  of  the  dynamic  stability  events  observed  during 
these  tests. 

The  ability  to  reliably  reproduce  these  grouped  extreme  waves  in  the  experimental  basin 
at  NSWCCD  will  enable  reduced  test  time  and  improved  repeatability  for  seakeeping  and 
dynamic  stability  tests.  This  technique  can  also  be  used  to  generate  improved  spectra  for 
irregular  wave  testing.  Because  of  the  temporal  and  spatial  limitations  of  the  basin,  it  is  difficult 
to  reproduce  the  entire  frequency  range  for  a  given  spectrum.  Generation  of  grouped  extreme 
waves  in  an  experimental  basin  will  enable  the  realistic  representation  of  more  severe  wave 
environmental  conditions  in  model  experiments,  which  is  impractical  with  current  irregular  wave 
testing  techniques. 
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Figure  2.  Wave  time-history  from  Chesapeake  Bay  Buoy  2,  on  10-20-2006  at  1414  hrs, 
significant  wave  height  as  reported  for  the  previous  30  min.  period. 
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Figure  3.  Wave  time-history  from  Chesapeake  Bay  Buoy  2,  on  10-20-2006  at  1445  hrs, 
significant  wave  height  as  reported  for  the  previous  30  min.  period. 
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Figure  4:  Wave  time-history  from  Chesapeake  Bay  Buoy  2,  on  10-23-2006  at  1030  hrs, 
significant  wave  height  as  reported  for  the  previous  30  min.  period. 


Figure  5:  Wave  time-history  from  Chesapeake  Bay  Buoy  2,  on  10-23-2006  at  1100  hrs, 
significant  wave  height  as  reported  for  the  previous  30  min.  period. 
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Theory 


Wave  Groups  with  Regular  Wave  Superposition 

Groups  of  asymmetric  extreme  waves  have  been  observed  in  nature,  but  have  not  been 
produced  experimentally.  For  this  investigation,  two  approaches  to  generate  extreme  waves  were 
studied:  the  wave-packet  method  and  finite-wave  linear  superposition.  The  wave-packet  method 
is  well  known  and  has  previously  been  demonstrated  as  a  method  to  produce  repeatable  single 
large-amplitude  waves,  at  a  fixed  location  in  an  experimental  basin.  To  generate  groups  of 
extreme  waves,  linear  superposition  of  a  series  of  finite  regular  waves  were  examined. 
Superposition  can  be  utilized  to  combine  wave  trains  of  variable  periods,  amplitudes,  and  cycles 
into  larger  or  smaller  amplitude  waves,  by  either  constructive  or  destructive  interference.  This 
technique  will  be  referred  to  as  the  finite-wave  linear  superposition  method. 

IVavc^Packel  Method 

The  underlying  principles  of  the  wave-packet  method  are  the  deep-water  dispersion 
relation  and  phase  velocity,  which  dictate  that  longer  waves  propagate  faster  than  shorter  ones. 
This  method,  previously  demonstrated  both  numerically  and  experimentally,  may  he  utilized  to 
generate  single  large  amplitude  waves.  At  the  upstream  position,  a  wave-packet  consisting  of  a 
band  of  frequencies  is  generated  by  a  wave-maker.  High  frequency  components  with  small 
amplitudes  are  generated  first  and  lower  frequency  components  arc  generated  subsequently.  As 
wave  components  propagate  downstream,  lower  frequency  components  coalesce  with  the  higher 
frequency  components  and  increasingly  larger  waves  are  formed.  In  general,  larger  waves  can  be 
produced  with  increased  distance  between  the  wave-maker  and  concentration  position. 

When  wave  components  are  superimposed  in-phase,  single  large  waves  or  symmetric 
groups  can  be  produced  in  irregular  seas.  The  characteristics  of  large  waves  arc  governed  by 
nonlinear  theory.  However,  linear  theory  was  used  for  the  first  stage  of  development.  A  wave- 
packet  ean  be  expressed  by  linear  superposition  of  the  wave  components, 

a(.xj)  =  Re^a„exp[/2^/„(/-/i()-/A„(.r-  v,)]  (I) 

n=() 


where  .v,  and  ta  represent  the  concentration  position  and  time,  respectively.  At  the  concentration 

position  and  time,  wave  components  are  superimposed  in-phase  and  a  large-amplitude  wave  may 
be  produced.  The  »th  frequency,  /„,  is  given  by 


(2) 


/.= 


AN 


where  A  represents  the  uniform  sampling  interval,  in  seconds,  and  N  is  the  total  number  of  discrete 

signals,  which  was  assumed  to  be  an  even  number  for  simplification.  The  duration  of  the  signal,  /,  is 
given  by  t  =  AN  and  the  signals  are  given  at  discrete  times  tn  =  nA  for  n  =  0,  1,  N  -  \ . 

The  wave  number,  kn ,  is  related  to  the  frequency,  /„,  by  the  deep  water  dispersion  relation, 


k  = 


 (2 njy 
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(3) 


For  a  wave-packet  that  propagates  in  the  downstream  (positive  x)  direction,  negative 
frequencies,  corresponding  to  -N/2<n<-\,  are  not  included.  Equation  (1)  describes  the 
wave-field  of  a  wave-packet,  at  a  given  position  and  time.  A  Rayleigh  distribution  can  be 
employed  as  a  normalized  spectrum, 

f2 

./  n 


-  f„ 

sn  =  —exp 
<7~ 


2(7 2 


(4) 


Figure  6  shows  the  spectrum  amplitude,  sn  f  computed  from  equation  (4),  for  rr  =  0.2  ^  N  =  8 1 92 

and  A  =  1/40  seconds.  Since  the  Nyquist  critical  frequency,  f  ,  is  20  Hz  (fNn  =  fc  = — - — ), 

(2A) 

the  frequency  range  is  from  0  to  20  Hz. 
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Figure  6.  Spectrum,  sn  when  cr  =  0.2 
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To  produce  large  waves  with  the  wave-packet  method,  the  important  parameters  are 
concentration  position,  .v, ,  concentration  time,  ta ,  and  spectrum,  sn.  The  concentration  position 

is  limited  by  the  length  of  basin,  the  concentration  time  may  be  arbitrarily  chosen  for  a  fixed 
concentration  position,  and  the  shape  and  width  of  the  amplitude  spectrum  may  be  selected  to 
provide  sufficient  energy  in  the  relevant  frequency  range.  Given  these  three  parameters,  the 
wave-field  is  described  by  equation  (1),  and  the  amplitude,  steepness,  and  breaking  point  of  the 
large-amplitude  wave  can  also  be  estimated.  Different  types  of  large-amplitude  waves  may  be 
obtained  with  combinations  of  different  amplitude  spectra  and  concentration  positions. 

Figure  7  shows  a  numerical  prediction  of  the  amplitudes  for  a  wave-packet  at  various 
stream-wise  positions  when  the  concentration  position  and  time  are  x  =  60  m  and  /tJ  =  100  sec, 

respectively.  From  the  amplitude  spectrum  given  in  equation  (4).  the  amplitudes  were  calculated 
from  equation  (1),  at  positions  .v  =  0,  20,  40,  50,  60.  and  70  m.  High  frequency  components  were 
followed  by  lower  frequency  components  in  the  upstream  locations.  All  wrave  components  were 
superimposed  in-phase  and  a  single  large-amplitude  wave  was  formed  at  x  =  60  m,  the  desired 
concentration  position.  Because  long  (low  frequency)  waves  propagate  faster  than  short  (high 
frequency)  waves,  due  to  the  phase  velocity  relationship  (e  =  A/T),  the  long  wave  components 
coalesce  with  the  short  waves  at  the  concentration  position.  Then,  for  the  same  reason,  the  higher 
frequency  components  lag  behind  the  lower  ones  after  the  concentration  position.  At  various 
positions,  although  magnitudes  of  the  fast  Fourier  transform  (FFT)  of  the  time  series  are  the 
same,  the  phases  are  different. 

Figure  8  presents  numerical  predictions  of  the  wave-packet  amplitudes  at  times  /  =70, 
80,  90,  100,  1 10,  and  120  seconds  so  the  examine  wave  steepness  could  be  examined.  At  the 
concentration  time  the  steepness  of  the  wave  is  large.  Onee  an  amplitude  sealing  factor  is  chosen 
to  seale  the  normalized  amplitude  to  a  desired  dimensional  value,  the  wave-breaking  point  may 
be  estimated  by  examining  the  steepness  of  the  waves. 
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Figure  7.  Time-history  of  packet  amplitude  for  a  notional  numerical  example, 
at  x  =  0,  20,  40,  50,  60,  and  70  m  when  xa  =  60  m  and  ta  =  1 00  sec 
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Figure  8.  Time-history  of  wave  height  for  a  notional  numerical  example, 
at  t  =70,  80,  90,  100,  1 10,  and  120  sec  when  xa  =  60  m  and  ta  - 100 sec 
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Finite-  IVave  Unear  Superposition  Method 

The  finite-wave  linear  superposition  method  employs  the  interaction  between  sequential 
finite  length  regular  waves,  of  varied  amplitude  and  frequency,  to  superimposed  at  a  desired 
location  in  the  basin,  resulting  in  a  group  of  asymmetric  large-amplitude  waves.  For  a  given 
finite  wave-train,  consideration  of  the  wave  amplitude,  period,  and  number  of  cycles  enabled  the 
different  finite  regular  waves  to  superpose  as  a  group  of  extreme  waves  at  a  desired,  repeatable 
location  in  the  basin. 

This  theory  was  examined  numerically  by  superposition  of  four  sine  waves  of  varying 
amplitudes  and  frequencies:  x,  =  l.lsin(q),  x,  =  sin(0.9q),  x,  =  1.2sin(l.h,),  .v4  =sin(0.8q).  The 

four  sine  waves  were  combined  to  form  periodic  asymmetric  groups  of  three  large  waves 
(Figure  9). 


Figure  9.  Combined  signal  output  from  sine  wave  superposition,  two  sine  waves  combined: 
jc,  +  .v2  (top)  and  x,  +  x4  (middle),  all  four  waves  ( x,  +  x,  +  x,  +  x4 )  combined  to  form 

groups  of  three  large  waves  (bottom) 
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In  this  example,  four  sine  waves  of  infinite  length  were  superposed,  and  the  three-wave 
group  was  shown  to  occur  periodically.  However,  four  finite-wave  sequences  would  result  in  a 
three-wave  group  occurring  at  only  one  location,  which  could  be  varied  based  on  the  individual 
finite-wave  amplitudes  and  periods. 

Linear  wave  theory  was  applied  to  calculate  the  group  velocity  and  regular  wave  groups 
were  calculated  to  coalesce  at  a  determined  concentration  position.  The  sequence  tiles  for  wave 
generation  were  determined  using  equations  (5-10). 

T 

Start  time  for  /7th  regular  wave  train:  tn  =t0+  —  ~  Tn (CVc  -  1 )  (5) 

4 


Voltage  for  /7th  wave  train  at  time  t :  V(t)=  H n  cos 


(6) 


Time  for  the  /7th  regular  wave  to  reach  test  point:  t r  =  — 1 —  +  — 

C,„  4 

Deep-water  group  velocity  of  /7th  regular  wave:  c .  =  <^~  = 

'  <Jk  An 

Wavelength,  A  =  T  2 
2  n 


Wave  steepness,  s  = 


// 

T 


(7) 

(8) 
(0) 
(10) 


where  T„  is  the  period  of  the  /7th  regular  wave  in  seconds,  //„  is  the  voltage  of  the  /7th  regular 
wave  (amplitude  of  flap  motion),  Cyc„  is  the  number  of  cycles  of  the  /7th  regular  wave,  .v„  is  the 
distance  to  the  concentration  location  (in  meters),  and  to  is  the  /ero  time  before  waves  start  in 
seconds.  Local  g  is  9.801 00  m/s2  ±0.0004. 

Wave  Croups  in  an  Irregular  Seaway 

Because  the  emphasis  of  this  study  was  the  generation  of  asymmetric  groups  of  extreme 
waves,  the  wave-packet  method  was  not  considered  for  this  part  of  the  study.  Previous  work  has 
shown  the  wave-packet  technique  to  be  successful  for  generating  single  large-amplitude  waves, 
or  symmetric  wave  groups,  and  embedding  single  extreme  waves  into  irregular  seas.  Much  of 
the  previous  work  has  focused  on  single  large-amplitude  waves  and  the  ability  to  generate  steep 
asymmetric  wave  groups  was  not  investigated.  For  this  investigation,  the  finite-wave  linear 
superposition  method  was  developed  and  employed  to  generate  asymmetric  groups  of  extreme 
waves  and  embed  them  into  two  scaled  irregular  sea  spectra,  a  Brctschncidcr  sea  state  8  and 
Hurricane  Camille. 
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Experimental  Methods  and  Test  Procedures 

Experiments  were  conducted  in  August  and  November  2007  in  the  Maneuvering  and 
Seakeeping  (MASK)  basin  at  NSWCCD  to  investigate  the  feasibility  of  producing  grouped 
extreme  waves  in  an  experimental  basin.  The  first  experiment  (Phase  I)  consisted  of  producing 
several  combinations  of  finite  regular  waves  with  varying  parameters,  including  amplitude, 
frequency,  and  signal  duration.  A  second  experiment  (Phase  II)  was  conducted  and  consisted  of 
two  parts.  The  first  part  applied  the  wave-packet  method  to  produce  single  large-amplitude 
waves,  based  on  the  technique  developed  by  Clauss  and  Bergmann  (1986)  and  Clauss  and 
Kuehnlein  (1994,  1995).  The  second  part  employed  the  achievements  from  Phase  1  to  embed 
grouped  extreme  waves,  obtained  through  finite  regular  wave  superposition  in  two  scaled 
irregular  seaways:  a  Bretschneider  sea  state  8  spectrum  and  a  Hurricane  Camille  spectrum. 
Wave-maker  settings  used  in  Phase  1  and  Phase  11  testing  are  detailed  in  Appendix  A. 
Descriptions  of  the  MASK  basin,  wake-maker  operations  for  regular  and  irregular  waves,  and 
two  wave  height  measurement  systems,  the  Senix  wave  gage  and  Global  Laser  Rangefinder 
Profilometry  (GLRP)  is  provided  in  the  following  sections. 

MASK  Description 

The  wave  generation  experiments  were  conducted  in  the  Maneuvering  and  Seakeeping 
(MASK)  basin  at  NSWCCD  (Figure  10).  Eight  pneumatic  wave-maker  units  are  located  along 
the  73  m  (240  ft)  side  of  the  basin  and  thirteen  units  along  the  110  m  (360  ft)  side  of  the  basin. 
The  basin  is  6  m  (20  ft)  deep.  A  1  15  m  (376  ft)  bridge  traverses  the  basin  and  can  be  moved  to  a 
45  degree  offset  from  the  longitudinal  center  of  the  basin.  The  two  perpendicular  banks  of  wave- 
makers  can  be  operated  individually  to  produce  long-crested  waves,  or  simultaneously  to 
generate  a  bi-directional  wave-field.  Sloping,  perforated,  concrete  beaches  are  located  on  each  of 
the  sides  of  the  basin  opposite  the  wave-makers  to  minimize  wave  reflections. 

Two  methods  are  used  to  control  the  flow  of  energy  into  the  wave-field:  varying  blower 
motor  speeds  supplying  air  to  the  pneumatic  domes  and  motion  amplitude  variation  of  the 
flapper  valve  that  controls  air  being  pumped  in  and  out  of  the  domes.  Hydraulic  cylinders  with  a 
±  10V  control  signal  are  employed  to  actuate  the  flapper  valves.  The  wave-maker  also  has  a 
series  of  lips  on  each  of  the  pneumatic  domes,  which  can  be  set  in  a  position  of  either  up  or 
down.  The  lips  can  be  used  to  modify  high  frequency  disturbances  in  the  generated  wave-field. 

For  these  experiments,  the  MASK  bridge  was  located  in  the  middle  of  the  basin,  parallel 
to  the  long  bank.  Wave  data  collected  from  sonic  probes  on  the  bridge  were  recorded  at  20  Hz, 
after  being  filtered  with  fixed  10  Hz  low-pass  filters.  Zeroes  were  taken  at  the  beginning  of  each 
testing  session  to  obtain  more  accurate  test  data  and  to  account  for  small  changes  in  the  water 
level  of  the  basin. 
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Wave-maker  Operation 

To  generate  long  crested  waves,  individual  wave-maker  units  in  each  bank  are  operated 
in-phase  and  produce  wave  segments  of  the  same  nominal  amplitude.  Wave  periods  can  be  either 
regular  or  random.  For  this  test,  waves  were  generated  from  the  short  bank  wave-makers  on  the 
left  side  and  travel  from  West  to  East  across  the  basin,  as  shown  in  Figure  10. 
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Regulag  Waves 

To  generate  regular  waves  the  wave-maker  control  program  employs  a  sequence  file  of 
voltages  programmed  at  40  cycles  per  second.  The  computer  generated  digital  control  voltage 
sequences  are  created  with  a  digital-to-analog  (D/A)  converter.  The  digitized  voltage  signal  is 
filtered  with  a  low-pass  smoothing  filter  and  a  high-pass  cutoff  frequency  prior  to  being  input 
into  the  wave-maker. 

In  this  study,  a  sequence  file  is  a  summation  of  four  finite  regular  waves  converted  to  a 
voltage  signal  for  the  hydraulic  wave-maker  piston.  Each  of  the  four  regular  waves  was 
controlled  by  blower  rpm,  maximum  voltage  (the  amplitude  of  flapper  motion),  frequency,  and 
the  number  of  wave  cycles.  Other  variables  for  the  calculation  of  a  wave  sequence  were  the  zero 
time  at  start,  which  indicates  a  series  of  zero  volts  before  initializing  the  wave,  and  the 
concentration  point,  which  is  the  point  at  which  the  waves  will  coalesce  to  form  a  desired  group 
of  three  large-amplitude  waves. 

Wave  sequences,  varying  four  separate  regular  wave  cycles  for  particular  pneumatic 
wave-maker  blower  rpm,  signal  amplitudes,  and  frequencies  (Table  1)  were  applied  to  generate 
wave  trains  designed  for  the  desired  groupiness,  with  repeatable  wave  properties,  at  a  fixed 
location  within  the  MASK.  The  finalized  sequence  file  was  a  summation  of  the  voltages  for  the 
four  wave  trains.  In  some  sequences,  a  flat  spot  existed,  due  to  the  signal  frequency  interaction, 
and  was  removed  manually  to  smooth  the  motion  of  the  hydraulic  pistons. 


Tabic  1.  Experimental  Parameters 


Parameter 

Value 

Blower  (rpm) 
Amplitude  (V) 
Frequency  (Hz) 
Number  of  cycles 

1100,  1300,  1500 

±  1.0,  2.0,  6.0,  8.0,  9.0,  9.5 

0.4-3. 3 

1-5 
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Irregular  Waves 

A  computer  program  was  used  to  generate  irregular  waves.  The  program  produced  digital 
control  signal  sequences  and  supplied  them  to  the  wave-maker  controller  through  a  digital-to- 
analog  (D/A)  converter.  The  digital  control  sequence  software  uses  a  filtered  white  noise 
technique.  To  produce  a  desired  wave  spectrum,  a  computer  generated  white  noise  signal  is  used 
to  actuate  the  flapper  valves,  resulting  in  wave-maker  dome  pressure  fluctuations.  The  wave 
energy  distribution,  as  a  function  of  the  frequency,  is  adjusted  with  the  driving  frequencies  for 
the  valves  controls.  The  blower  speed  is  adjusted  to  control  the  wave  amplitude.  To  represent 
nominal  conditions  for  extreme  waves,  two  long-crested  irregular  wave  spectra  were  considered 
for  this  study,  a  Bretschneider  wave  spectrum  for  sea  state  8  and  a  Hurricane  Camille  wave 
spectrum.  Both  speetra  were  modeled  using  two  different  scale  ratios.  Plots  of  the  theoretical 
Brctsehneider  spectral  energy  distribution  for  sea  state  8  and  the  Hurricane  Camille  spectral 
energy  distribution,  which  was  actually  measured  during  the  hurricane,  arc  shown  in 
Appendix  C. 

Senix  Wave  Cages 

To  measure  the  generated  waves,  the  MASK  basin  contains  an  array  of  non-contact 
ultrasonic  sensors,  suspended  from  the  MASK  Bridge  (Table  2).  Locations  (.vv)  of  the  sensors  are 
referenced  from  the  south-east  comer  of  the  MASK  basin  (Figure  10).  The  Senix  Corporation 
model  ULTRA-SR-BP  sonic  sensor  transmits  an  ultrasonic  wave  and  measures  the  time  of 
reflection  from  the  target  to  calculate  the  distance.  Wave  height  measurements  are  collected  from 
an  array  of  six  bridge-mounted  wave  probes  (Figure  10).  The  probes  arc  connected  to  a  computer 
on  the  bridge  that  transmits  the  wave  data  to  a  shore  computer,  where  the  data  arc  collected  and 
stored. 

Each  probe  (Figure  1 1 )  emits  a  conical  sonic  beam  with  a  nominal  12-dcgrce  total  angle. 
Measurements  for  steep  waves  may  have  drop-outs  because  the  water  surface  tends  to  scatter  the 
sonic  beam  away  from  the  sensor.  When  attempting  to  measure  large,  steep  w'aves,  the  probes 
can  have  signal  drop-outs  or  inaccurate  measurements  due  to  scattering,  which  may  result  in 
measurement  error.  Calibration  results  and  uncertainty  estimates  for  the  tests  arc  described  in  the 
next  section  of  the  report. 
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able  2.  Wave  Height  Sonics  Located  on  MASK  Bridge 
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Figure  11.  Senix  ULTRA-SR-BP  wave  gage  and  attachment  on  the  MASK  bridge 

GLRP 

Global  Laser  Rangefinder  Profilometry  (GLRP)  was  developed  at  NSWCCD  to  provide  a 
eapability  for  a  time-resolved  field  measurement  of  wave  elevations.  GLRP  provides  high-rate, 
three-dimensional  mapping  of  surfaee  waves  over  a  large  physical  area.  It  is  suitable  for  the 
study  of  large-amplitude  transient  wave  dynamics  in  experimental  facilities.  The  GLRP  concept 
and  original  prototypes  were  demonstrated  and  documented  previously,  and  a  fully-funetional 
GLRP  system  was  constructed  in  the  Maneuvering  and  Seakeeping  (MASK)  basin  at  NSWCCD 
(Atsavapranee,  et  al.,  2005;  Cameal,  et  al.,  2005a). 

The  MASK  GLRP  system  was  designed  to  provide  aeeurate  model-seale  whole-field 
measurements  of  wave  surfaee  profiles  during  free-running  and  captive  model  tests.  The  GLRP 
system  illuminates  the  water  surfaee  with  distinct  points  using  laser  diodes  and  measures  vertical 
fluctuations  (Figure  12).  The  apparent  position  of  each  diode  was  calibrated  and  recorded  on  a 
charge-eoupled  device  (CCD)  detector.  The  eapability  of  the  system  to  measure  waves  created  in 
the  MASK  facility,  as  well  as  the  wave-field  generated  by  a  sealed  ship  model,  was 
demonstrated  in  previous  testing  (Carneal,  et  al.  2005b;  Cameal  and  Atsavapranee,  2006).  The 
absolute  uncertainty  of  the  system  was  estimated  at  3  mm,  whieh  represents  a  relative 
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uncertainty  on  the  order  of  0.5%  of  the  full  measurement  range  (Cameal,  et  al.,  2005b;  Cameal 
and  Atsavapranee,  2006).  Calibration  results  and  uncertainty  estimates  for  the  tests  are  described 
in  the  next  section  of  the  report. 

The  GLRP  system  implemented  in  the  MASK  facility  is  shown  in  Figure  13.  The  system 
used  for  this  experiment  consisted  of  two  panels,  with  200  diodes  in  each  panel,  spanning  an  area 
of  1.5  m  by  3  m  (4.92  ft  by  9.84  ft),  at  a  spatial  resolution  of  7.62  cm  (19.35  in.)  in  each 
direction.  Two  CCD  cameras,  with  a  frequency  of  30  Hz  and  a  resolution  of  1392  X  1040  pixels, 
are  also  mounted  on  each  panel.  Each  panel  is  mounted  to  the  GLRP  platform  with  an  ACME 
screw  and  two  linear  bearings  to  allow  traversal  of  the  entire  panel.  The  platform  has  a  seeding 
system  consisting  of  two  seeding  rakes,  with  three  nozzles  each,  to  disperse  fluorescent  dye 
across  the  measurement  area.  The  platform  contains  all  the  necessary  hardware  for  seeding  and 
data  acquisition  (Figure  14).  A  more  complete  description  of  the  system,  calibration  procedures, 
and  data  processing  methods  used  in  the  MASK  GLRP  system  is  contained  in  Cameal,  et  al. 
(2005)  and  Cameal  and  Atsavapranee  (2006). 


Figure  12.  GLRP  array  illuminating  water  surface  in  MASK,  as  viewed  from  the  MASK 

Bridge 
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Figure  13.  GLRP  array  panel  in  MASK 
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Instrumentation  Calibration  and  Uncertainty 


The  following  section  describes  the  calibration  of  the  Senix  Ultra-SR-BP  wave  gages  and 
Global  Laser  Rangefinder  Profilometry  (GLRP)  and  the  uncertainty  estimates  of  the  calibrations. 

Senix  Ultrasonic  Wave  Gages 

The  Senix  wave  gages  were  calibrated  in-situ,  by  varying  the  distance  above  a  measured 
calm  water  level  in  the  MASK.  The  transducers  were  located  via  precision  machined  pin 
locations  with  an  estimated  uncertainty  of  ±0.127  mm  (±0.005  inch).  The  calibration  range  was 
±381  mm  (±15  inches),  in  127-mm  (5-inch)  increments.  The  calibration  locations  were  -381,  - 
254,  -127,  0,  ±127,  ±254,  ±381  mm  (-15,  -10,  -5,  0,  ±5,  ±10,  ±15  inches). 

The  analog  output  for  the  Senix  is  0  to  10  Vde.  Data  from  the  Senix  were  collected  and 
processed  on  a  PC.  The  data  were  digitized  with  a  National  Instruments  16-bit  data  acquisition 
card,  N1  6036,  at  a  sample  rate  of  24  Hz.  The  anti-aliasing  filters  are  Frequency  Devices 
D78L8L,  set  at  a  cutoff  frequency  of  8  Hz.  Typical  data  collection  time  was  at  least  60  s. 

\ormal  ( 'a!  i  brut  ion 

The  normal  nominal  gain  setting  for  the  Senix  is  76.2  mm/V  (3  mehes/V).  The  noise 
level  of  the  Senix  for  measurements  over  calm  water  is  typically  between  1  and  10  mV.  The 
maximum  Type  A  expanded  uncertainty  from  the  ISO  GUM  (1995)  is  ±0.048  mm  (±0.0019 
inch)  for  1,000  samples,  or  41.7  s  averaging  time.  The  Type  A  standard  uncertainty  is  defined  as 

uA=^  (11) 

where  s/,  is  the  standard  deviation  of  the  height  and  n  is  the  number  of  samples.  The  expanded 
uncertainty  is  defined  as 


U  =  kfu  (12) 

where  kf  is  the  coverage  factor.  Normally  at  the  95%  confidence  level,  kf=  1.959  ~  2  .  The 
maximum  combined  and  expanded  uncertainty  from  these  two  elements  is  0.135  mm  (0.0053 
inch).  The  contribution  of  the  Senix  noise  is  negligible,  compared  to  the  machining  accuracy  of 
±0.127  mm  (±0.005  inch)  for  the  positioning  device. 

An  example  calibration  plot  is  presented  in  Figure  15,  as  a  residual  plot,  for  Senix  wave 
gage  #1 .  A  residual  plot  is  the  difference  between  the  data  and  a  linear  regression  fit.  The  dashed 
line  in  Figure  15  is  the  prediction  limit  at  the  95%  confidence  level,  from  the  calibration  theory 
of  Scheffe  (1973)  and  Carroll,  et  al.  (1988).  The  uncertainty  in  the  measurement  from  the 
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instrument  noise  and  position  is  smaller  than  the  symbols  in  the  figure.  The  red  symbol  in  the 
figure  is  an  outlier  and  was  excluded  from  the  linear  regression  analysis. 

The  outlier  is  caused  by  the  over-ranging  of  the  Senix,  in  this  case  at  10  V.  Due  to  the 
location  of  the  Senix  sensor  relative  to  the  water  at  0  (zero)  mm,  the  transducer  may  over-range 
at  either  0  or  10  V,  but  not  both.  In  a  few  cases,  the  sensor  was  located  relative  to  the  water  so 
that  over-ranging  did  not  oeeur. 

Most  of  the  calibration  uncertainty  is  due  to  the  data  scatter  in  the  calibration.  The 
contributions  from  instrument  noise  and  position  are  quite  small  by  comparison.  The  relative 
uncertainty  in  the  calibration  is  presented  in  Figure  16.  As  Figure  16  indicates,  the  maximum 
relative  uncertainty  is  about  ±2%.  The  relative  uncertainty  is  computed  from  the  full-range  (fr) 
value,  which  is  nominally  381  mm  ( 15  inches). 

Table  3  is  a  summary  of  the  calibration  data  acquired  before  the  test.  In  this  ease,  the 
calibration  data  were  reported  for  the  position  of  the  probe  above  the  water.  In  the  experiment, 
the  gage  remained  fixed,  and  the  water  moved  relative  to  the  gage.  In  the  data  processing  routine, 
the  signs  of  the  gain  and  slope  were  opposite  those  in  the  table.  The  slope  is  negative  and  the 
intercept  positive  for  the  wave  measurements. 

In  Table  3,  the  column  SEE  is  the  standard  error  estimate  from  the  linear  regression 
analysis,  while  the  r  in  the  next  column  is  the  correlation  coefficient  The  SEE  is  essentially  the 
standard  deviation  of  the  curve  tit  and  is  a  measure  of  the  randomness  of  the  data.  The  precise 
definition  may  be  found  in  most  texts  on  statistical  analysis,  sueh  as  Ross  (2004).  For  a  highly 
linear  eurve  fit,  r  is  near  one  as  indicated  in  the  table,  and  SEE  will  be  near  zero. 

The  column  t-test,  is  the  result  of  a  hypothesis  test  for  a  comparison  of  the  measured 
slope  to  the  nominal  slope  of  76.2  mm/V  (3  inehes/V)  from  Ross  (2004).  The  column  %  is  the 
value  of  the  inverse  Student-/  from  the  hypothesis  test.  Any  value  from  the  t-test  that  is  less  than 
tt/s  is  statistically  the  same  as  the  nominal  slope.  From  the  table,  the  only  wave  gage  that  is 
different  from  the  nominal  slope  is  #6. 

The  last  two  columns  are  the  expanded  uncertainty  of  the  gages  at  the  95%  confidence 
level,  at  10  Vde.  The  last  eolumn  is  percent  full-range  (%  fr)  from  the  measured  slope,  where  the 
nominal  full-scale  value  is  381  mm  (15  inches).  As  the  table  indicates,  the  largest  uncertainty 
was  about  ±8  mm  (±0.3 1  inches),  or  ±2%  fr,  at  the  95%  confidence  level. 
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Higher  Range  Calibration 

To  measure  larger  waves,  the  gain  for  the  Senix  wave  gages  was  ehanged  to  a  nominal 
value  of  101.6  mm/V  (4  inches/V).  However,  the  gages  were  not  calibrated  until  after  the  Phase 
II  test.  The  results  described  here  are  from  the  post-test  calibration.  The  calibration  was  limited 
by  the  calibration  fixture,  over  a  total  range  of  762  mm  (30  inches).  The  required  range  for  0  to 
10  Vdc  was  ±508  mm  (±20  inches).  From  the  limitations  of  the  fixture,  the  7  calibration 
locations  were  then  -508,  -381,  -254,  -127,  0,  ±127,  and  ±254  (-20,  -15,  -10,  -5.  0,  +5,  and  +10 
inches). 

The  calibration  results  for  Senix  wave  gage  #1  are  presented  in  Figure  17.  The  measured 
uncertainty  for  each  point  is  again  smaller  than  the  symbols.  The  noise  level  for  the  Senix  was 
0.8  to  6.2  mV,  which  was  similar  to  the  noise  level  from  the  normal  range  calibration.  The  Type 
A  uncertainty  had  a  maximum  expanded  uncertainty  of  ±0.040  mm  (±0.0016  inch).  The  resulting 
combined  uncertainty  was  then  ±0.13  mm  (±0.0052  inch).  The  uncertainty  from  the  instrument 
noise  was  small  in  comparison  to  uncertainty  of  the  probe  location,  and  most  of  the  uncertainty 
was  in  the  curve  fit. 

As  Figure  17  indicates,  the  maximum  uncertainty  was  about  ±17  mm.  The  uncertainty  at 
+508  mm  (+20  inches)  was  significantly  higher  than  at  -508  mm  (-20  inches).  The  larger  value 
is  from  the  extrapolation  of  the  calibration  data  from  +254  mm  (+10  inches)  to  +508  inches  (+20 
inches).  The  relative  uncertainty  is  shown  in  Figure  18.  The  maximum  relative  uncertainty  for 
this  gage  was  ±3.4%  fr. 

The  calibration  results  are  summarized  in  Table  4.  The  data  trends  are  similar  to  the 
normal  calibration  (Tabic  3).  All  of  the  probes,  with  the  exception  of  #6,  have  a  slope  that  is 
statistically  the  same  as  the  nominal  gain  setting  of  101.6  mm/V  (4  inches/V).  With  the 
exception  of  gage  #1,  the  relative  uncertainty  of  the  gages  was  less  than  or  equal  to  ±2%  fr, 
similar  to  previous  results.  The  cause  of  the  larger  uncertainty  in  probe  #1  is  unknown. 
Calibration  over  the  full  range,  ±508  mm  (±20  inches),  of  the  Senix  wave  gages  may  reduce  the 
uncertainty. 

Because  both  the  normal  and  higher  range  calibrations  are  performed  on  a  fiat-plane 
surface,  the  aeeuracy  associated  with  signal  drop-outs  when  attempting  to  measure  steep  waves 
should  be  considered.  However,  the  influence  of  wave  steepness  on  signal  drop-outs  for  the 
Senix  probes  is  not  currently  known. 
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u95  (%)  Wave  Height  Residuals  (mm) 
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Figure  15.  Calibration  of  Senix  Wave  Gage  #1  for  Normal  Range 
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Fable  3.  Summary  of  Senix  Wave  Gage  Calibrations  for  Normal  Range 


■ 

Probe # 

Slope 

Intercept 

SEE 

r 

t-test 

^95 

U95 

u95 

(mm  A/) 

(mm) 

(mm) 

(mm) 

(%  fr) 

1 

76  138 

-392.67 

1.410 

0  999979 

0.308 

245 

7  33 

1.93 

3 

76  325 

-380  34 

1.808 

0.999975 

0.611 

2.36 

7  99 

2  09 

4 

76  269 

-390  20 

0.945 

0.999991 

0.515 

245 

4.92 

1.29 

5 

76  314 

-381.16 

1.844 

0.999977 

0.545 

2  45 

8  50 

2.23 

6 

76  538 

-392  73 

0.679 

0.999995 

3.51 

2.45 

3.53 

0  92  ! 

8 

76.313 

-395.16 

0.358 

0.999999 

2.24 

2.45 

1  85 

0  49  1 

Table  4.  Summary  of  Senix 

Wave  Gage 

Calibrations  for  High  Range 

Probe  # 

Slope 

Intercept 

SEE 

r 

t-test 

t95 

u95 

U95 

(mm/V) 

(mm) 

(mm) 

(mm) 

(%fr) 

1 

102  477 

-512  02 

2  821 

0  999956 

2  06 

2.57 

17.29 

3  38 

3 

101.458 

-491  95 

1.054 

0.999992 

0.704 

2.78 

7.53 

1.48 

4 

101  442 

-512.87 

1.221 

0  999992 

0.857 

2.57 

7.33 

1.44 

5 

101  883 

-502.24 

1.469 

0.999985 

1.002 

2.78 

10.47 

2.06 

6 

102.260 

-511.20 

0.621 

0.999997 

5.56 

2.78 

4.38 

0  86 

8 

101.635 

-515.52 

0.622 

0.999997 

0.290 

2  78 

4.35 

0  86 

GLRP 

Calibration  of  the  GLRP  measurement  was  performed  by  raising  and  lowering  the  GLRP 
panels  to  several  standoff  heights.  Images  of  the  laser  diode’s  interseetion  with  the  water  surface 
were  recorded  at  eaeh  height  to  simulate  vertical  motion  of  the  water  surface.  This  process  is 
termed  surfaee  height  calibration,  and  a  superimposed  calibration  image  for  a  single  eamera  is 
shown  in  Figure  19.  In  this  experiment,  calibration  heights  ranged  from  -304.8  mm  to  +304.8 
mm  (-1  ft.  to  +12  ft.)  in  increments  of  25.4  mm  (1  in.).  The  calibration  data  provided  the 
positions  of  the  images  of  the  laser  diodes  on  the  CCD  eamera  in  pixels  versus  the  height  or  z- 
coordinate.  A  typieal  result  is  shown  in  Figure  20.  Results  are  presented  for  two  resolutions  of 
the  measurements.  Calibration  data  for  this  report  were  based  upon  7  points  (Figure  20a),  at  an 
increment  of  101.6  mm  (4  in.).  The  full  data  set  is  shown  in  Figure  20b,  at  a  resolution  of  25.4 
mm  ( 1  in.). 

The  calibration  for  this  report  was  obtained  from  a  seeond-order  polynomial  fit.  Figure 
20a  indicates  that  a  seeond-order  fit  may  not  be  justified.  A  measure  of  the  uncertainty  is  the 
Standard  Error  of  Estimate  (SEE).  As  the  figure  indicates,  the  SEE  is  significantly  redueed  with 
an  increased  number  of  data  points.  For  7  points,  the  SEE  was  nearly  19  mm,  compared  to 
12.1  mm  for  23  points,  where  one  point  was  excluded  as  an  outlier  from  the  24  point  data  set. 
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The  data  were  then  evaluated  for  an  optimal  curve  fit  with  a  reasonably  lower  value  of 
SEE.  The  results  are  summarized  in  Table  5  for  diode  #80,  which  had  one  of  the  lower  values  of 
SEE  in  the  data  set  of  80  diodes.  As  the  table  indicates,  a  third-order  polynomial  fit  has  a 
significantly  lower  value  of  SEE,  4.6  mm,  compared  to  an  SEE  of  10  mm  for  the  second-order 
fit,  with  7  points  for  the  data  in  this  report. 

Table  5.  Comparison  of  Polynomial  Fits  of  Different  Orders  and  Number  of  Points  for 
GLRP  Diode  #80  Calibration  for  Wave  Height 


No.  Points 

Order 

SEE 

(mm) 

r 

7 

1 

18.6 

0.9970 

7 

2 

10.4 

0.9993 

23 

1 

12.1 

0.9981 

23 

2 

7.8 

0.9992 

23 

3 

4.6 

0.9997 

The  calibration  result  for  diode  #80,  as  a  third-order  polynomial,  is  presented  in  Figure  21 
as  a  residual  plot.  As  the  figure  indicates,  on  the  basis  of  the  prediction  limit  at  the  95% 
confidence  level,  the  uncertainty  was  10  mm,  which  is  similar  to  the  high  range  calibration  of 
Senix  gage  #5  in  Table  4.  Senix  gage  #5  was  the  closest  to  the  GLRP.  However,  the  uncertainty 
varies  with  diode  number,  or  location,  as  indicated  by  the  SEE  in  Figure  22.  The  SEE  decreased 
systematically  from  a  value  of  14.6  mm  at  diode  #  1  to  4.6  mm  at  diode  #80.  The  value  from  the 
prediction  is  about  2.2  *SEE.  Thus,  the  uncertainty  in  the  wave  height  calibration  varied  between 
10  mm  and  32  mm.  The  variation  of  the  3rd-order  polynomial  coefficients  with  diode  number  is 
shown  in  Figure  23.  These  coefficients  also  vary  systematically  with  diode  number,  with  trends 
similar  to  SEE. 

In  addition  to  calibration  of  the  system  in  the  r-coordinate,  an  .vy-mapping  calibration  was 
also  performed.  Because  the  camera  was  at  an  angle  to  the  water  surface,  the  xy-positions  of  the 
diodes  must  be  determined  as  a  function  of  water  height.  This  was  accomplished  by  recording 
images  of  a  precision  machined  calibration  target  at  several  heights,  at  a  known  position  in  the 
MASK  facility.  A  sample  image  is  shown  in  Figure  24.  These  images  were  then  de-warped  using 
standard  image  processing  techniques.  The  calibration  positions  of  each  diode  image  were 
calculated  from  a  corresponding  x-  and  ^-position  in  the  MASK  facility,  for  each  diode  at  each 
height.  A  second-order  polynomial  fit  forx  andy,  as  a  function  of  distance  from  neutral  position, 
was  computed  for  each  measurement  point.  The  second-order  polynomial  fits  were  obtained  for 
x,  y,  z  physical  coordinates  in  the  MASK  facility,  as  a  function  of  diode  position  on  the  CCD 
detector. 
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Data  processing  was  accomplished  by  blob  analysis  of  the  raw  data  images  and  locating 
the  diode  image  in  each  frame.  Information  from  the  calibration  process  was  used  to  assign  the 
detected  blobs  to  the  appropriate  calibration  curves.  The  centroid  of  the  detected  blob  was 
determined  and  the  calibration  curves  were  applied  to  determine  the  location  (x,  y,  z)  of  each 
blob  in  each  frame,  and  the  data  were  stored.  Because  the  eurrent  MASK  system  generates  a 
large  amount  of  raw  data,  approximately  170MB/s  in  a  four  eamera,  two  panel  configuration,  the 
data  processing  must  be  performed  offline.  A  more  complete  description  of  the  system, 
calibration  procedures,  and  data  processing  methods  in  the  MASK  GL.RP  system  is  documented 
in  Atsavapranee,  et  al.  (2005),  Cameal,  et  al.  (2005),  and  Cameal  and  Atsavapranee  (2006). 


Figure  19.  Superimposed  surface  height  calibration  image 
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Figure  20.  GLRP  Linear  Wave  Height  Calibration  for  Diode  #80 
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Figure  24.  Sample  image  for  XY  calibration 
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Experimental  Results 


Wave  Croups  with  Regular  Wave  Superposition 

The  wave-paeket  method  and  finite-wave  linear  superposition  method  were  used  to 
generate  extreme  waves.  The  wave-paeket  method  is  well  known  and  has  been  demonstrated  as  a 
method  to  produee  repeatable  single  large-amplitude  waves,  at  a  fixed  location  in  an 
experimental  basin.  Earlier  research  did  not  investigate  or  indicate  the  capability  of  this  method 
to  generate  realistic  groups  of  asymmetric  extreme  waves.  To  develop  this  later  capability,  the 
finite-wave  linear  superposition  method  was  used.  Through  linear  superposition  of  these  finite 
regular  waves,  unidirectional  asymmetrie  groups  of  extreme  waves  were  produced  at  repeatable, 
fixed  locations  in  the  basin.  Regular  wave  results  from  the  wave-paeket  method  and  the  finite- 
wave  linear  superposition  method  are  shown  in  Appendix  B.  Results  from  embedded  wave 
groups  in  irregular  waves  are  shown  in  Appendix  C.  GLRP  results  and  comparisons  to  Senix 
Probe  measurements  are  shown  in  Appendix  D. 

lVu}j>Packcl  Method 

The  wave-paeket  method  applies  wave  components,  which  arc  superimposed  in-phase,  to 
produee  single  large  waves  in  irregular  seas.  T  his  method  was  applied  to  generate  a  single  large- 
amplitude  wave  at  a  fixed  location  in  the  MASK  basin.  Linear  superposition  of  components 
using  Rayleigh  distributions  as  a  speetrum,  coupled  with  varied  scale  parameters  and  amplitudes, 
were  used  to  generate  large-amplitude  waves  for  this  part  of  the  study.  For  comparison, 
additional  measurements  were  taken  of  pure  sine  waves  generated  with  varied  frequencies  and 
amplitudes.  The  test  results  with  a  single  frequency  sinusoidal  input  signal  and  with  a  wave- 
paeket  input  signal  are  diseussed.  Results  are  shown  in  Appendix  B. 

As  a  baseline  test,  a  sine  wave  with  the  frequency  of  0.3  Hz  was  given  as  an  input  for  the 
wave-maker.  Figure  25  shows  the  time  series  data  from  the  six  wave  gages  installed  along  the 
bridge  in  the  basin.  The  positions  of  the  wave  gages  are  shown  in  Figure  10.  In  general,  the 
measured  frequencies  of  the  Senix  signals  were  in  good  agreement  with  the  input  frequency. 
However,  the  amplitude  of  the  sine  wave  measured  at  probe  8,  located  at  the  most  upstream 
position,  was  much  smaller  than  the  other  downstream  probe  signals.  The  wave  amplitude  at 
probe  4  was  a  little  larger  than  at  probe  3.  The  generated  wave-field  in  the  model  basin  is  not 
completely  uniform  (O’Dea  and  Newman,  2007;  Smith,  et  al.,  2007). 

Linear  theory  indicates  that  if  an  upstream  wave-paeket,  for  example  at  x  =  0,  can  be 
generated,  then  a  large  wave  at  a  desired  eoneentration  position,  at  v  =  60  m,  ean  be  expeeted 
(Figure  7).  The  input  signal  for  the  wave-maker  ean  be  determined  if  the  transfer  function  of  the 
wave  maker  is  known.  However,  the  transfer  function  of  the  MASK  wave-maker  is  not  clearly 
known,  particularly  the  phase.  Additionally,  in  many  cases,  the  transfer  functions  are  not  linear. 
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Because  the  wave-maker  transfer  function  of  the  MASK  is  not  well  known,  the  signal  at 
x  =  0  (Figure  7)  was  utilized  as  an  input  to  the  wave-maker.  The  probe  data  are  plotted  in  Figures 
25  and  26.  The  probe  signals  indicate  that  the  concentration  position  was  located  between  probes 
3  and  1 .  The  maximum  wave  height  measured  at  probe  3  is  about  1 8  cm  and  at  the  upstream 
probe  8  is  7.65  cm.  Additional  measurements  between  probes  3  and  1  are  necessary  to  determine 
the  exact  location  of  the  concentration  position  and  the  maximum  wave  height. 

To  evaluate  the  effect  of  the  amplitude  on  the  wave-maker  signal,  three  input  signals  with 
different  amplitudes  were  supplied  to  the  wave-maker.  The  maximum  values  of  the  wave-maker 
flapper  control  signals  were  3,  6,  and  9  V,  shown  in  Figure  27  by  the  blue,  green,  and  red  lines, 
respectively.  The  measured  probe  data  are  plotted  in  Figure  28.  In  these  figures,  the  probe  data 
were  multiplied  by  scale  factors  of  the  input  signals,  3  for  blue,  1.5  for  green,  and  1  for  red.  The 
probe  data  from  three  different  input  amplitudes  agree  quite  well. 

The  results  (Figure  28)  demonstrate  that  linear  theory  can  identify  an  optimum  input 
signal.  Various  input  signals,  with  different  amplitude  spectra,  can  produce  a  desired  large  wave 
signal,  at  a  given  position.  Once  an  input  signal  is  selected,  the  signal  can  be  amplified  to 
produce  a  larger  wave,  as  long  as  a  breaking  wave  does  not  occur  in  the  wave-packet  during  the 

process.  Linear  theory  also  shows  that  a  wave-field  at  any  position  and  time  can  be  determined 
by  a  few  parameters:  concentration  position,  xu ,  concentration  time,  ta ,  and  amplitude  spectrum, 

an .  Because  concentration  position  is  limited  by  the  length  of  the  model  basin  and  concentration 

time  may  be  chosen  arbitrarily,  the  amplitude  spectrum  is  the  only  free  parameter  Comparisons 
of  the  amplitudes  for  different  voltages  and  the  resulting  waves  are  shown  in  Figures  29  and  30. 

Results  of  this  study  indicate  that  a  database  with  different  wave-maker  input  signals  is 
an  essential  first  step  in  producing  single  large-amplitude  waves  deterministically.  This  database 
will  identify  an  input  signal  that  will  produce  a  desired  wave-packet,  at  a  specified  position.  The 
shape  of  the  wave-packet  should  be  the  primary  focus.  The  test  results  show  that  a  larger  wave 
can  be  produced  by  multiplying  an  input  signal  by  a  desired  scale  factor.  Because  the  generated 
wave-field  in  the  MASK  was  shown  to  be  repeatable,  large  waves  may  be  produced  with  the 
identified  input  signals.  The  results  of  this  study  confirm  that  a  wave-field  may  be  modeled  well 
with  linear  theory. 

Finite-  Wave  Linear  Superposition  Method 

Based  on  the  principle  of  superposition,  regular  waves  of  finite  length  and  varying 
amplitude  were  considered  as  a  method  to  produce  groups  of  asymmetric  large-amplitude  waves. 
For  this  part  of  the  study,  fifty-six  runs  were  completed  in  the  MASK,  The  influence  of  the 
different  parameters  were  examined  by  thirty-six  different  wave  sequence  files  and  a  repeatable 
group  of  large-amplitude  waves  were  obtained  at  a  fixed  location  in  the  basin.  Time-histories  of 
the  waves,  recorded  by  the  Senix  probes,  are  shown  in  Appendix  C. 


38 


A  number  of  runs  containing  large  wave  groups  were  observed.  Some  sample  runs:  15, 
35,  43,  51,  and  61  are  presented  where  the  wave  sequences  produced  coalesced  to  form  a  group 
of  three  waves  of  successively  increasing  amplitudes  that  occur  at  probes  3  and  4  (Figure  3 1 ).  In 
some  cases,  the  third  wave  in  the  group  reached  a  wave  height  of  nearly  0.76  m  (30  inches) 
(Figure  32).  The  wave  characteristics  of  these  groups,  obtained  from  the  wave-time  histories 
recorded  by  the  sonic  probes,  arc  shown  (Figures  41-55).  Some  of  the  observed  trends  are 
discussed  below. 

The  wave  time-histories  show  a  wave-packet  generated  by  the  wave-maker,  as  defined  in 
Table  1 1,  arriving  at  probe  8.  The  varying  speeds  of  the  waves  within  the  packet  caused  the 
waves  to  form  the  desired  three-wave  group  formation  within  the  region  measured  by  probes  4 
and  3,  before  dissipating  toward  the  end  of  the  basin,  as  shown  in  the  measurements  from 
probe  1  (Figures  41-55). 

From  the  conditions  tested,  an  increase  in  the  blower  rpm  led  to  an  overall  increase  in  the 
generated  wave  amplitudes.  Setting  the  lips  on  the  wave-maker  to  the  down  position  resulted  in  a 
slight  improvement  in  the  overall  shape  of  the  waves,  because  the  lips  filter  out  high  frequency 
disturbances.  For  the  range  of  parameters  tested,  a  generalized  limit  for  the  wave-maker  blower 
rpm  was  established.  For  wave  cycles  with  large  amplitudes,  indicated  by  large  voltage  signals,  a 
blower  setting  of  1 300  rpm  could  not  be  exceeded  without  breaking  waves  occurring.  For  some 
cases  (ex.  Run  61),  reduced  signal  voltages  for  the  1st  and  4th  cycles  enabled  a  blower  setting  of 
1500  rpm  to  be  applied  without  breaking  waves  resulting. 

As  observ  ed  in  runs  15  and  43,  a  high  voltage  signal  was  applied  for  the  first  of  the  wave 
cycles.  This  was  combined  with  small  amplitude  signals  and  a  low  number  of  cycles  for  the 
second  and  third  waves.  For  the  fourth  wave  cycle  signal,  there  was  a  decrease  in  the  amplitude 
and  increase  in  the  frequency.  The  combination  of  these  four  waves  resulted  in  the  occurrence  of 
a  larger  and  steeper  third  wave  in  the  group,  at  probe  4  in  the  basin. 

As  observed  in  runs  15  and  35,  a  small  increase  in  the  third  and  fourth  wave  cycle 
frequencies,  with  other  parameters  constant,  resulted  in  a  larger  amplitude,  but  slightly  less  steep 
second  wave,  at  probe  4.  As  observed  in  runs  35  and  43,  a  small  increase  in  frequency  and 
decrease  in  amplitude  of  the  third  wave  cycle,  and  a  small  decrease  in  the  frequency  and  increase 
in  amplitude  of  the  fourth  wave  cycle,  resulted  in  larger  amplitude  and  steeper  second  waves  at 
probe  4.  However,  neglecting  changes  to  the  third  and  fourth  wave  cycles  most  likely  results  in  a 
larger  and  steeper  third  wave.  Because  the  wave  characteristics  exceeded  the  calibration  limits  of 
the  Senix  probes  and  the  steepness  of  the  largest  wave,  this  was  unconfirmed  during  this  set  of 
experiments  due  to  saturation  and  drop-outs  in  the  sonic  probe  measurements. 

As  observed  in  runs  51  and  61,  a  decrease  in  the  first  and  fourth  cycle  amplitudes,  with 
one  additional  cycle  added  to  the  second  wave  signal,  resulted  in  a  greater  wave  height  for  the 
third  wave  in  the  group,  still  with  1/10  wave  steepness.  This  decrease  also  resulted  in  larger  and 
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steeper  first  and  second  waves  in  the  wave  group. 

At  probe  3,  many  of  the  runs  showed  a  decrease  in  the  first  and  second  wave  heights  of 
the  three-wave  group  compared  to  the  wave  heights  of  the  group  measured  at  probe  4.  For  the 
third  wave  in  the  group,  the  largest  wave  amplitude  was  typically  observed  at  probe  4.  However, 
the  greatest  steepness  of  the  third  wave  in  the  group  was  typically  measured  at  probe  3.  At 
probe  3,  the  second  wave  in  the  group  was  slightly  less  steep,  but  the  third  wave  was  slightly 
steeper  than  the  same  group  at  probe  4.  This  difference  in  the  wave  group  characteristics  is 
important  because  it  provides  variations  in  wave  heights  and  steepness  at  two  locations  close 
together  in  the  basin,  to  enable  a  model  test  to  accurately  simulate  a  range  of  conditions. 

As  part  of  Phase  II  testing,  runs  5-8  and  42,  were  performed  to  repeat  the  optimal  wave 
groups  obtained  through  the  finite-wave  linear  superposition  method  in  Phase  1.  The  wave  time- 
histories  of  the  wave  groups  were  in  good  agreement  between  the  August  and  November  2007 
tests  (Tables  6  and  7).  Recalibration  and  gain  change  of  the  Senix  probes  enabled  measurements 
of  larger  wave  amplitudes  without  saturation,  although  some  runs  still  had  wave  amplitudes 
larger  than  ±  50  cm,  which  remained  beyond  the  recalibrated  range  of  the  Senix  probes. 


Table  6.  Summary  of  Maximum  Wave  Height  Results  for  Finite-Wave  Linear 


Superposition  Method-  Probe  #3 


Run  # 

h(m) 

h/A 

h/L 

hfS  (m)** 

15 

0.615 

0.125 

0.157 

28.7 

7* 

0.738 

0.143 

0.189 

34.4 

35 

0.595 

0.111 

0.152 

27.7 

42* 

0.678 

0.111 

0.173 

31.6 

43 

0.606 

0.125 

0.155 

28.2 

8* 

0.744 

0.143 

0.190 

34.7 

51 

0.633 

0.111 

0  162 

29  5 

5* 

0.633 

0.111 

0.162 

29.5 

61 

0.595 

0.125 

0.152 

27.7 

6* 

0.705 

0.143 

0.180 

32.9 

*Run  numbers  from  Phase  II  repeat  runs 

**  For  a  scale  ratio  of  46.6 
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Table  7.  Summary  of  Maximum  Wave  Height  Results  for  Finite-Wave  Linear 


Superposition  Method-  Probe  #4 


Run  # 

h(m) 

h/A 

h/L 

hfs  (m)** 

15 

0.743 

0  100 

0.190 

34  6 

7* 

0.739 

0.091 

0  189 

34.4 

35 

0.722 

0  100 

0.185 

33  6 

42* 

0.669 

0.077 

0.0  171 

31.2 

43 

0  706 

0.091 

0  180 

32.9 

8* 

0.732 

0.091 

0  187 

34.1 

51 

0  666 

0.100 

0  170 

31.0 

5* 

0.507 

0.100 

0  130 

23.6 

61 

0  727 

0.111 

0.185 

33.9 

6* 

0.672 

0.100 

0.172 

31.3 

*Run  numbers  from  Phase  11  repeat  runs 
**  For  a  scale  ratio  of  46.6 

Wave  Groups  in  an  Irregular  Seaway 

These  experiments  were  performed  as  part  of  Phase  II,  and  consisted  of  extending  the 
methodology  for  generating  large-amplitude  wave  groups  using  the  finite-wave  linear 
superposition  method,  developed  in  Phase  I.  In  the  Phase  II  experiments,  the  large-amplitude 
wave  groups  were  embedded  in  two  different  irregular  seaways.  A  list  of  the  conditions  tested  is 
shown  in  Table  12.  Throughout  this  phase,  additional  measurements  of  the  wave-field  were 
acquired  with  the  Global  Laser  Rangefinder  Profilometery  (GLRP)  system.  Results  are  shown  in 
Appendix  C. 

Scnix  Results 

Two  wave  spectra  were  examined  to  determine  the  feasibility  of  the  embedding  method, 
a  30th  scale  Bretschneider  Sea  State  8,  with  Hs=  38.1  cm  (15.0  in.)  and  Tm=  3.0  s,  and  two 
Hurricane  Camille  spectra,  a  30th  scale,  with  Hs=  40.64  cm  (16.0  in.)  and  T„,=  2.45  s,  and  a  46th 
scale,  with  Hs=  26.16  cm  (10.30  in.)  and  Tm=  1.96  s.  The  two  scales  considered  in  this  study 
represent  typical  scale  ratios  for  ship  model  testing  at  NSWCCD,  used  to  assess  seakeeping  and 
dynamic  stability  performance. 

As  a  first  comparison,  run  15  from  the  Phase  1  testing  was  embedded  at  both  the  end  and 
in  the  middle  of  the  30th  scale  Bretschneider  Sea  State  8  spectra.  In  both  cases  the  wave  group 
retained  the  general  characteristics  after  being  embedded.  Because  of  the  desired  application  to 
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dynamic  stability  model  testing,  the  process  was  continued  with  the  wave  group  embedded  in  the 
middle  of  the  irregular  seaway  time  series  (Figures  61-71).  Results  from  Senix  probe 
measurements  at  probes  3  and  4  for  the  embedded  large-amplitude  wave  group  are  shown  in 
Tables  X  and  9. 

To  assess  the  impact  of  the  embedded  wave  group  on  the  spectral  shape,  spectral  analysis 
was  performed  of  runs  with  both  the  Bretschneider  and  Hurricane  Camille  spectra 
(Figures  94-96).  Examination  of  the  spectra  shows  a  good  agreement  between  the  idealized  and 
generated  spectra  with  the  embedded  wave  group  for  the  30th  scale  Bretschneider  Sea  State  X. 
With  the  cmdedded  wave  group,  the  peak  of  the  spectra  was  significantly  higher,  but  the  overall 
spectral  shape  was  still  close  to  the  idealized  value.  However,  the  idealized  and  generated  spectra 
for  the  two  scaled  Hurricane  Camille  seaways  did  not  agree  as  well.  The  30th  scale  spectrum 
with  the  embedded  wave  groups  had  higher  spectral  values  at  the  lower  frequencies  and  a  lower 
peak  value  compared  to  the  idealized  spectrum.  The  46th  scale  Hurricane  Camille  spectrum  with 
the  embedded  wave  group  had  a  higher  peak  than  the  idealized  spectrum,  but  also  a  secondary 
peak,  near  a  frequency  of  3.X  rad/s.  Because  a  deterministic  time  series  was  used  for  the 
Hurricane  Camille  sequence,  it  is  possible  that  not  enough  of  the  spectral  components  were 
included  to  define  a  realistic  spectral  shape,  skewing  the  peak  of  the  spectral  shape.  With  the 
extreme  wave  group  embedded  in  the  time  series,  a  secondary,  larger  peak  may  have  occurred, 
more  closely  approximating  the  peak  of  the  idealized  Camille  spectrum.  Additional  testing  will 
have  to  be  performed  to  confirm  this  conjecture. 

Variation  of  the  embedded  wave  groups  in  irregular  seaways  produced  runs  where  three- 
wave  groups  occurred  in  the  seaway  with  two  outcomes.  In  the  first,  successively  larger  waves 
led  to  a  steep  wave,  which  could  be  classified  as  an  extreme  wave  (H  >  2.0  Hs).  Iri  the  second,  a 
wave  of  larger  amplitude  than  the  previous  two  in  the  group,  and  at  least  two  waves  still 
exceeding  the  criteria  for  an  extreme  wave,  were  observed.  Both  of  these  outcomes  were 
reproducible  and  repeatable  (Figures  61-71).  Processing  of  the  characteristics  of  the  wave 
groups,  as  each  of  the  three  waves  occurred  in  the  irregular  seaway,  was  completed 
(Figures  72-93). 
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Table  8.  Summary  of  Maximum  Wave  Height  Results  for  Finite-Wave  Linear 
Superposition  Method  Embedded  Into  Irregular  Waves-  Probe  #3 


Run  # 

Spectrum* 

h(m) 

h/A 

h/L 

hfS  (mj** 

ii 

BS  SS8 

0  804 

0.143 

0.205 

37.5 

14 

BSSS8 

0.799 

0.143 

0.204 

37.2 

21 

BS  SS8 

0.602 

0.111 

0.153 

28.1 

22 

BS  SS8 

0.644 

0.111 

0.164 

30  0 

24 

BS  SS8 

0  718 

0.125 

0.184 

334 

26 

HC 

0  647 

0.018 

0.165 

30.2 

27 

HC 

0.530 

0.071 

0.136 

24.7 

30 

HC 

0.542 

0.015 

0.139 

25.3 

41 

HC 

0.480 

0.063 

0.123 

22.4 

43 

HC 

0.649 

0.091 

0.165 

30  2 

Table  9.  Summary  of  Maximum  Wave  Height  Results  for  Finite-Wave  Linear 


Run  # 

Spectrum* 

h  (m) 

h/A 

h/L 

hfs  (m)** 

11 

BS  SS8 

0.714 

0.091 

0.183 

33.3 

14 

BS  SS8 

0.710 

0.100 

0.182 

33.1 

21 

BS  SS8 

0.525 

0.063 

0  134 

24.5 

22 

BS  SS8 

0.530 

0.050 

0.136 

24.7 

24 

BS  SS8 

0.540 

0.526 

0.138 

25.2 

26 

HC 

0.761 

0.040 

0.195 

35.5 

27 

HC 

0.582 

0.067 

0.149 

27.1 

30 

HC 

0.786 

0.100 

0.201 

36.6 

41 

HC 

0.478 

0.040 

0.122 

22.3 

43 

HC 

0.494 

0.048 

0.126 

23.0 

*Embedded  into  Bretschneider  sea  state  8  (BS  SS8)  and  Hurricane  Camille  (HC)  spectra 

**  For  a  scale  ratio  of  46.6 
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GLRP  Results 


Results  from  GLRP  measurements  and  comparisons  to  Senix  probe  data  are  shown  in 
Appendix  D.  Three-dimensional  representations  of  GLRP  data  are  presented  every  third  figure, 
in  Figures  97-121.  These  representations  show  several  snapshots  of  the  wave  group  evolution  in 
the  measurement  region.  The  first  frame  shows  the  water  surface  before  the  wave  group  enters 
the  measurement  region.  The  second  frame  shows  the  beginning  of  the  wave  group.  In  the  third 
and  fourth  frames,  the  maximum  and  minimum  wave  heights  encountered  in  the  measurement 
are  shown.  The  final  frame  shows  the  calm  water  surface  at  the  end  of  the  run. 

A  single  point  time-history  for  a  selected  GLRP  diode  is  shown  every  third  figure,  in 
Figures  98-122.  The  time-history  clearly  shows  the  evolution  of  the  wave  group,  from  the 
occurrence  of  the  smaller  waves  at  the  beginning  to  the  larger  waves  in  the  middle  of  the  wave 
group.  At  the  end  of  the  measurement,  the  water  surface  returns  to  a  relatively  calm  state  after 
the  wave  group  has  traversed  the  measurement  area.  In  addition,  sonic  measurements  were  taken 
simultaneously  with  GLRP  measurements.  The  sonic  probe  in  closest  proximity  to  the  GLRP 
platform  was  probe  5.  A  comparison  between  the  time  histories  of  a  single  GLRP  point  and  sonic 
probe  5  is  shown  every  third  figure,  in  Figures  99-123.  Qualitatively  the  measurements  agree 
very  well.  Because  the  systems  were  not  spatially  co-located  in  the  basin,  a  truly  quantitative 
comparison  cannot  be  made.  However,  the  amplitudes  appear  to  be  within  the  uncertainty  in  the 
measurements,  and  most  of  the  difference  may  be  the  phase  shift  due  to  location.  These  details 
have  not  been  evaluated  for  this  data,  in  future  tests,  a  sonic  probe  should  be  placed  within  the 
measurement  region  of  the  GLRP  platform  to  provide  a  direct  comparison  between  the  two 
measurement  techniques. 

During  Run  7  (Figures  97-99),  the  GLRP  system  began  collecting  data  after  the  smallest 
waves  in  the  wave  group  had  already  reached  the  measurement  area.  The  evolution  of  the  wave 
group  is  apparent,  as  the  interaction  between  waves  produced  minimum  and  maximum  wave 
heights  occurring  at  approximately  1 7  and  1 8  seconds,  respectively.  After  the  maximum  wave 
arrived,  the  wave  height  fluctuations  decayed  quickly,  and  the  basin  began  to  settle.  The  GLRP 
measurement  showed  a  similar  peak,  with  a  slightly  lower  minimum  than  the  sonic 
measurements,  as  well  as  a  faster  decrease  in  wave  activity  after  the  maximum  wave  encounter 
occurred. 

In  Run  9  (Figures  103-105),  a  selected  portion  of  the  time  series  is  shown  to  emphasize 
the  evolution  of  the  beginning  of  the  wave  group.  For  this  case,  a  set  of  grouped  extreme  waves 
were  added  to  the  wave-field  after  the  development  of  Bretschneider  sea  state  8.  A  comparison 
between  the  time  histories  of  a  single  GLRP  point  and  sonic  probe  5  is  shown  (Figure  105)  and 
the  GLRP  and  sonic  measurements  show  a  strong  qualitative  agreement. 

For  Run  1  1,  and  all  successive  runs,  data  are  shown  for  Cameras  1  and  2  only.  The  data 
from  Camera  3  exhibited  a  high  noise  level.  Data  could  be  extracted  from  Camera  3  with  further 
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software  modification,  but  it  was  decided  to  replace  the  camera  for  future  tests  instead  of 
exhausting  resources  to  develop  unnecessary  high  noise  level  data  processing  applications. 

A  single  point  time-history  for  a  selected  GLRP  diode  and  Senix  probe  5  is  shown 
(Figure  108).  For  this  case,  waves  were  generated  to  simulate  a  Bretschneidcr  sea  state  8  with  a 
large-amplitude  wave  group  embedded  in  the  wave-field.  The  minimum  and  maximum  wave 
heights  occurred  at  approximately  52  and  54  seconds.  The  GLRP  and  sonic  measurements  show 
a  strong  qualitative  agreement.  The  GLRP  system  measured  a  lower  trough  and  peak  than  the 
sonic  probe,  but  these  differences  could  stem  from  the  difference  in  measurement  locations  and 
capabilities  of  the  two  systems. 

A  single  point  time-history  for  a  selected  GLRP  diode  is  shown  in  Figure  1 16,  which 
shows  a  section  of  the  time-history  in  the  middle  of  the  wave  generation  sequence.  For  this  ease, 
waves  were  generated  to  simulate  conditions  during  Hurrieane  Camille  with  a  wave  group 
embedded  into  the  wave-field.  The  maximum  and  minimum  wave  heights  occurred  at 
approximately  61  and  62  seconds  after  the  GLRP  system  data  collection  began.  The  large- 
amplitude  wave  group  is  clearly  evident  in  the  middle  of  the  GLRP  time  series,  surrounded  by 
waves  occurring  in  the  simulated  hurrieane  conditions.  A  comparison  between  the  time-histories 
of  a  single  GLRP  point  and  sonic  probe  5  is  shown  in  Figure  117.  In  this  ease,  the  maximum 
wave  height  encounters  measured  by  the  GLRP  and  sonic  systems  were  very  similar,  while  the 
GLRP  system  continued  to  deteet  a  deeper  trough  for  the  largest  wave  than  the  sonic  probe. 

A  single  point  time-history  for  a  selected  GLRP  diode  is  shown  in  Figure  122.  For  this 
ease,  waves  were  generated  to  simulate  the  conditions  during  Hurricane  Camille  with  a  wave 
group  embedded  in  the  wave-field.  The  maximum  wave  encounter  occurred  at  approximately  57 
seconds  after  GLRP  system  data  collection  began.  A  comparison  between  the  time  histories  of  a 
single  GLRP  point  and  sonic  probe  5  is  shown  in  Figure  123.  The  GLRP  and  sonic 
measurements  show  a  strong  qualitative  agreement.  As  in  previous  measurements,  the  GLRP 
system  measured  a  lower  trough  and  peak  than  the  sonic  probe.  Only  a  few  peaks  and  troughs 
exhibited  major  differences  in  this  run. 

The  GLRP  measurements  performed  in  this  work  were  very  successful  in  obtaining  3-D 
spatially  dense  measurements  of  the  wave-fields.  Comparisons  to  Senix  probes  show  very  strong 
qualitative  agreement  between  the  measurement  techniques.  The  observed  differences  between 
the  measurement  techniques  can  be  attributed  to  several  factors.  The  most  important  factor  is  that 
the  measurements  were  taken  in  different  locations  in  the  basin,  since  the  sonic  probe  is 
approximately  4.57  m  (15  ft)  closer  to  the  short  bank  wave-maker  than  the  closest  GLRP  point. 
Given  that  the  wave-fields  generated  in  this  experiment  ranged  from  sinusoidal  to  highly 
random,  the  qualitative  agreement  between  the  GLRP  and  sonic  measurements  was  expected  to 
decrease  with  additional  random  wave  interaction.  The  data  collected  in  this  experiment 
supported  this  expectation,  and  for  the  isolated  large-amplitude  wave  groups  (runs  7  and  8),  the 
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qualitative  agreement  between  the  measurements  was  very  favorable.  For  the  Bretschneider  sea 
state  8  and  Hurricane  Camille  test  runs,  the  qualitative  agreement  between  the  measurement 
systems  was  less  favorable,  especially  in  regions  of  random  wave  interaction  away  from  the 
wave  groups.  However,  the  measurements  display  a  level  of  qualitative  agreement  that  was  very 
close,  even  in  the  random  seaway  cases,  despite  the  different  location  from  probe  5. 

In  future  tests,  a  sonic  probe  will  be  placed  in  the  same  location  as  the  GLRP  system  to 
provide  a  quantitative  comparison  between  the  measurement  systems.  In  addition,  directional 
wave  probes  will  be  deployed  to  determine  the  ability  of  GLRP  measurements  to  provide 
directional  wave  spectra.  In  order  to  exactly  locate  the  sonic  probes,  laser  measurement  and 
alignment  systems  will  be  used  to  precisely  determine  the  sonic  locations  in  the  basin,  as  well  as 
the  exact  position  of  the  GLRP  platform  and  measurement  areas  as  a  function  of  grid  angle. 
Future  work  may  also  include  reprocessing  of  the  data  contained  in  this  report  to  extract 
directional  wave  information  from  the  collected  data  to  determine  the  capability  of  the  current 
GLRP  system  in  the  measurement  of  directional  wave  spectra.  Additionally,  the  differences  in 
amplitudes  of  the  waves  from  both  systems  should  be  compared  quantitatively  to  the  estimated 
uncertainties.  The  phase  shift  in  the  measurements  should  be  checked  for  consistency  in  the 
relative  location  of  the  two  types  of  sensors. 

Comparison  ofSenix  qnd_  GUjP  Results 

Results  including  the  measured  maximum  wave  heights  from  the  Senix  Probe  #5  and 
GLRP  point  measurements  for  probes  within  at  x=3.4-3.5  m,  y=  1.3- 1.4  m,  at  one  of  the  middle 
points  on  the  first  panel.  The  percent  difference  for  each  run  are  shown  in  Table  10.  The  percent 
difference  was  calculated  using  the  following  standard  formula: 


%  Diff  = 


I^Soin  K, 


GLRP] 


h  Senix  +  h(;LRP 


•100 


(13) 
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Tabic  10.  Summary  of  Comparison  of  Maximum  Wave  Height  Results  GLRP  and  Senix 


Probe  #5  Measurements 


Run# 

Senix  h 

(m) 

GLRP  h 

(m) 

Difference 

(%) 

7 

0.629 

0  651 

3.44 

8 

0.613 

0.634 

3.37 

11 

0.653 

0  635 

2.80 

21 

0  561 

0.573 

2  12 

22 

0.626 

0  629 

0.478 

24 

0.661 

0  616 

7  05 

26 

0.635 

0  618 

2.71 

27 

0.435 

0.482 

10.25 

30 

0.552 

0.494 

11.09 

All  of  the  data  in  Table  1 0  are  within  the  limits  of  the  Senix  calibration  (h=  1 .0 1 6  m),  but 
some  exeeed  the  calibration  range  of  the  GLRP  system  (H=  0.610  m).  Additional  work  is 
neeessary  to  determine  the  calibration  limit  of  GLRP.  Runs  7.  9,  11,21,  22,  and  26  in  Table  10 
are  within  the  measurement  uncertainty  of  Senix  probe  #5  and  probably  have  a  significant 
overlap  of  uncertainty  estimates  between  the  two  systems. 


Summon'  of  Results 

The  maximum  calibrated  attained  wave  height,  as  normalized  with  a  nominal 
representative  ship  model  length  (h/L),  was  0.205  compared  to  0. 1  for  the  typieal  largest  regular 
wave  height  to  ship  model  length  ratio  used  for  testing.  This  is  double  the  maximum  wave  height 
normally  achieved  in  the  MASK.  The  maximum  wave  steepness  (h/>.)  observed  was 
approximately  1/7,  approaching  the  theoretical  limit  for  a  non-breaking  wave,  compared  to  a 
1/10  wave  steepness  vvhieh  is  typieally  the  maximum  for  regular  wave  testing.  The  largest 
estimated  full-seale  wave  height  observed  was  37.5  m  (123  ft.)  at  a  46.6  seale  ratio. 
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Conclusions  and  Recommendations 


The  purpose  of  this  investigation  was  to  demonstrate  the  feasibility  of  producing  extreme 
waves,  both  singularly  and  in  groups,  in  an  experimental  basin  for  applications  to  surface-ship 
model  experiments.  Physical  model  experiments  provide  a  valuable  evaluation  tool  to  assess 
surface  ship  performance  in  extreme  conditions.  Currently,  the  best  representation  of  a  seaway  in 
the  model  basin  is  through  the  modeling  of  spectral  shapes  by  generating  random  waves.  Model 
experiments  conducted  in  random  waves  to  assess  dynamic  stability,  slamming,  and  ultimate 
strength  require  very  long  run  times  to  ensure  that  the  extreme  wave  events,  which  are  probable 
for  the  seaway,  have  been  encountered.  Employing  a  deterministic  model  testing  technique  will 
reduce  testing  time  in  the  model  basin,  enhance  the  realistic  representation  of  the  seaway,  and 
increase  confidence  in  the  assessment  of  dynamic  stability  and  secondary  structural  loads 
performance  in  extreme  conditions. 

Grouped  extreme  waves  of  varied,  yet  repeatable,  amplitudes  and  locations  were 
produced  during  the  course  of  this  investigation.  The  NSWCCD  MASK  basin  pneumatic  wave- 
maker  generated  extreme  wave  groups  in  the  basin.  Superposition  of  regular  waves  created  wave 
groups  and  embedded  them  in  an  irregular  seaway.  The  wave-packet  method  generated  single 
large-amplitude  waves.  The  generation  of  large-amplitude  wave  groups  was  demonstrated  in  the 
basin,  with  the  goal  of  eventually  producing  simulated  extreme  wave  groups  in  irregular  waves, 
similar  to  those  observed  in  nature.  GRLP  measurements  were  shown  to  provide  a  similar  level 
of  accuracy  as  measurements  using  the  Senix  wave  probes,  but  also  enabled  3D  wave-field 
measurements,  in  addition  to  point  measurements,  of  very  large,  steep  waves.  This  study 
demonstrated  the  initial  steps  for  the  development  of  a  sophisticated  and  efficient  experimental 
technique  for  dynamic  stability  model  testing  which  will  enable  the  realization  of  wave  groups 
within  a  given  seaway  that  will  induce  severe  motions. 

Future  Work 

The  work  outlined  in  this  report  details  progress  made  in  FY07  and  the  first  half  of  FY08. 
Expanding  on  the  initial  success  of  this  work,  the  eventual  objective  of  improved  dynamic 
stability  and  structural  testing  of  a  free-running  ship  model  with  wave  groups  in  the  MASK  may 
be  achieved.  Plans  for  future  work  are  presented  below. 

One  of  the  first  tasks  will  be  additional  investigation  of  the  Chesapeake  Bay  wave  buoy 
data.  The  wave  buoy  data  can  be  evaluated  to  develop  assessment  methods  to  determine  the 
frequency  of  occurrence  of  groups  of  extreme  waves,  better  characterize  these  wave  groups,  and 
possibly  correlate  the  groups  to  dynamic  stability  events.  A  criterion  to  identify  extreme  wave 
groups  could  be  applied  to  buoy  and  satellite  imagery  data  from  ocean  measurements,  for  several 
geographic  regions  and  time-periods.  This  would  aid  with  initial  development  to  determine  the 
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probability  of  oeeurrenee  for  groups  of  extreme  waves,  which  eould  be  expanded  to  assess  the 
global  oeean  wave  environment. 

A  significant  emphasis  of  additional  work  will  inelude  the  generation  of  bidirectional 
waves  from  the  short  bank  and  long  bank  in  the  MASK.  A  primary  goal  of  this  work  is  to  re¬ 
create  wave  time  histories  seen  in  nature,  such  as  those  measured  from  the  Chesapeake  Bay  wave 
buoys  or  other  oeean  measurements,  to  provide  realistie  extreme  wave  groups  for  application  in 
experimental  testing.  Investigation  of  techniques  for  embedding  additional  multiple  large- 
amplitude  waves  groups,  sueh  as  a  five  wave  group,  in  an  irregular  wave  spectrum,  in  both  sea 
states  7  and  8,  can  be  investigated.  Further  adjustment  of  wave  eyele  parameters,  to  aehieve  a 
range  of  grouped  wave  amplitudes  and  locations,  may  be  investigated.  Although  the  finite-wave 
linear  superposition  method  enabled  repeatable,  predictable  generation  of  three  wave  groups, 
further  application  of  the  wave-paeket  approach  should  also  be  examined  to  attempt  to  produee 
wave  groups  with  large-amplitude,  asymmetric,  steep  waves  and  embed  them  in  an  irregular 
seaway. 

The  application  of  the  methodology  for  dynamie  stability  testing  may  be  demonstrated  by 
testing  of  seale  models  with  the  wave  groups,  initially  for  dead-ship  conditions,  but  also  with 
forward  speed,  to  assess  dynamie  stability  performance.  In  addition  to  dynamie  stability 
performance  evaluation,  the  large-amplitude  wave  group  technique  can  be  utilized  to  evaluate 
secondary  loads  on  ships  and  offshore  structures.  A  theory  will  need  to  be  developed  for 
understanding  how  to  perform  statistical  extrapolation  of  ship  motion  response  from  extreme 
wave  groups  to  determine  lifetime  risks  for  ship  stability  failures,  perhaps  building  on  some  of 
the  reeent  work  of  Alford  (Alford,  et  al„  2006;  Alford.  2008). 

Continuing  throughout  this  future  work,  additional  measurements  of  wave-fields  will  be 
taken  with  the  Global  Laser  Rangefinder  Profilometery  (GLILP)  system.  Development  of  a 
Infrared  Global  Laser  Rangefinder  Profilometery  (IRGLRP)  system  which  does  not  require  dye 
seeding  of  the  water  surfaee  is  currently  underway.  Applications  of  this  system  to  large- 
amplitude  wave  group  testing  will  also  be  explored. 

The  calibration  range  of  the  measurement  systems  should  be  extended  to  match  the 
maximum  wave  heights  of  the  eurrent  research.  Some  of  the  wave  heights  in  the  experiments 
continued  to  exeeed  the  recalibrated  range  of  the  Senix  probes,  ±  0.508  m.  The  calibration  limit 
of  the  GLRP  system  has  not  yet  been  investigated.  Calibration  fixtures  for  both  systems  should 
be  modified  for  a  wave  height  measurement  range  larger  than  ±  0.6m 
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Appendix  A:  Wave-maker  Settings 


Tabic  1 1.  Phase  I  Wave  Run  Matrix 


Run 

Number* 

Cycle  1 

Cycle  2 

Cycle  3 

Cycle  4 

Blower 

RPM 

Lips 

\\  ave 

sequence  file 

02 

0  5  Hz 

0.4  Hz 

5  V 

5  V 

10  cycles 

1 0  cycles 

1300 

Down 

groupsOl  .t\t 

03 

2  Hz 

2.25  Hz 

2.5  Hz 

2.98  Hz 

9  5V 

9.5  V 

9.5  V 

9.5  V 

5  cycles 

1  cycle 

1  cycles 

1  cycle 

1300 

Down 

test9.dat 

04 

2  Hz 

2.25  Hz 

2.5  Hz 

2.98  Hz 

9  5  V 

9.5  V 

9.5  V 

9.5  V 

5  cycles 

1  cycle 

1  cycles 

1  cycle 

1300 

Down 

groups02.txt 

05 

2  Hz 

2.25  Hz 

2.5  Hz 

2.98  Hz 

9  5  V 

9.5  V 

9.5  V 

9  5  V 

5  cycles 

1  cycles 

1  cycles 

1  cycles 

1300 

Down 

groups03.txt 

06 

2  Hz 

2.25  Hz 

2.5  Hz 

2.98  Hz 

9  5  V 

7  V 

7  V 

9.5  V 

5  cycles 

1  cycles 

1  cycles 

1  cycles 

1300 

Down 

groups04.txt 

07 

2  Hz 

2.25  Hz 

2.5  Hz 

2.98  Hz 

9.5  V 

3  V 

3  V 

9.5  V 

5  cycles 

1  cycles 

1  cycles 

1  cycles 

1300 

Down 

groups05.txt 

09 

2  Hz 

2.25  Hz 

2.5  Hz 

2.98  Hz 

9.5  V 

1  V 

1  V 

9.5  V 

5  cycles 

1  cycles 

1  cycles 

1  cycles 

1300 

Down 

groups06.txt 

10 

2  Hz 

2.25  Hz 

2.5  Hz 

2.98  Hz 

9.5  V 

0.1  V 

0.1  V 

9.5  V 

5  cycles 

1  cycles 

1  cycles 

1  cycles 

1300 

Down 

groups07.txt 

1  1 

2  Hz 

2.25  Hz 

2.5  Hz 

2.98  Hz 

9  5  V 

0.1  V 

0  1  V 

9.5  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups08.txt 

13 

2  Hz 

2.25  Hz 

2.5  Hz. 

2.98  Hz 

9  5  V 

1  V 

1  V 

9.5  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups09.txt 
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14 

2  Hz 

6  V 

4  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

2.98  Hz 

4  V 

3  cycles 

1300 

Down 

groupsl0.txt 

15 

2  Hz 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

2.98  Hz 

9.5  V 

2  cycles 

1300 

Down 

groups09.txt 

16 

2  Hz 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

2.98  Hz 

9.5  V 

2  cycles 

1300 

up 

groups09.txt 

17 

2  Hz 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

2.98  Hz 

9.5  V 

2  cycles 

1300 

Up 

groups09.txt 

18 

2  Hz 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

2.98  Hz 

9.5  V 

2  cycles 

1000 

Up 

groups09.txt 

19 

2  Hz 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

2.98  Hz 

9.5  V 

2  cycles 

1 150 

Up 

groups09.txt 

20 

2  Hz 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

2.98  Hz 

9.5  V 

2  cycles 

1500 

Up 

groups09.txt 

21 

2.1  Hz 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

3.3  Hz 

9.0  V 

2  cycles 

1300 

Up 

groups  1 1  txt 

22 

2.1  Hz 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

3.3  Hz 

9.0  V 

2  cycles 

1100 

Up 

groups  1  l.txt 

23 

2.1  Hz 

9.5  V 

5  cycles 

2.25  Hz  1 

V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

3.3  Hz 

9.0  V 

2  cycles 

1100 

Up 

groupsl  l.txt 

24 

2.1  Hz 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

I  cycles 

2.5  Hz 

1  V 

1  cycles 

3.4  Hz  9.0 

V 

2  cycles 

1300 

Up 

groupsl2.txt 

25 

2.1  Hz 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

1  V 

1  cycles 

3.2  Hz 

9.0  V 

2  cycles 

1300 

Up 

groupsl3.txt 
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26 

2.1  H  7 

9.5  V 

5  cycles 

2.25  Hz  1 

V 

1  cycles 

2.6  Hz 

2  V 

1  cycles 

3.2  Hz  9.0 

V 

2  cycles 

1300 

Up 

groups  I4.txt 

27 

2.1  Hz 

2.25  Hz 

2.3  Hz 

3.2  Hz 

9.5  V 

1  V 

2  V 

9.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Up 

groups  1 5.txt 

29 

2.1  Hz 

2.25  Hz 

2.3  Hz 

3.2  Hz 

9.5  V 

1  V 

2  V 

9.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Up 

groupsl5.txt 

30 

2.1  Hz 

2.25  Hz 

2.3  Hz 

3.2  Hz 

9.5  V 

1  V 

2  V 

9.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1100 

Up 

groups  1 5.txt 

31 

2.1  Hz 

2.4  Hz 

2.3  Hz 

3.2  Hz 

9.5  V 

1  V 

2  V 

9.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1 300 

Up 

groups  16.txt 

32 

2.1  Hz 

2.25  Hz 

2.3  Hz 

3.2  Hz 

9.5  V 

2  V 

2  V 

9.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Up 

groups  1 7.txt 

33 

2.1  Hz 

2  25  Hz 

2.17  Hz 

3.2  Hz 

9.5  V 

2  V 

2  V 

9.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1 300 

LIT 

groupsl8.txt 

34 

2.1  Hz 

2.25  Hz 

2.3  Hz 

3.2  Hz 

9.5  V 

2  V 

2  V 

9.0  V 

5  cycles 

1  cycles 

2  cycles 

2  cycles 

1300 

Up 

groups  1 9.txt 

35 

2  Hz 

2.25  Hz 

2.7  Hz 

3  Hz 

9.5  V 

1  V 

1  V 

9.5  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Up 

groups20.txt 

36 

2  Hz 

2.25  Hz 

2.7  Hz 

3  Hz 

9.5  V 

1  V 

1  V 

9.5  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1500 

Up 

groups20.txt 

37 

2  Hz 

2.25  Hz 

2.7  Hz 

3  Hz 

9.5  V 

1  V 

1  V 

9.5  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1400 

Up 

groups20.txt 
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38 

2  H? 

9.5  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.7  Hz 

1  V 

1  cycles 

3.1  Hz 

9.5  V 

2  cycles 

1300 

Up 

groups2 1  .txt 

39 

2  H  7 

2.25  Hz 

2.5  Hz 

3.1  Hz 

9.5  V 

1  V 

1  V 

90  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Up  

groups22.txt 

40 

2  Hz 

2.25  Hz  1 

2.5  Hz 

3.1  Hz 

9.5  V 

V 

2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Up 

groups23.txt 

41 

2  H  7 

2.25  Hz 

2.5  Hz 

3.1  Hz 

9.5  V 

1  V 

2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups23.txt 

42 

2  H  7 

2.25  Hz 

2  5  Hz 

3.1  Hz 

9.5  V 

1  V 

2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups24.txt 

43 

2  H7 

2.25  Hz 

2.5  Hz 

3.1  Hz 

9.5  V 

1  V 

2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups25.txt 

2  H  7 

2.25  Hz 

2.5  Hz 

3.1  Hz 

44 

9.5  V 

1  V 

2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1 100 

Down 

groups25.txt 

2  H7 

2.25  Hz 

2.5  Hz 

3.1  Hz 

45 

9.5  V 

1  V 

2.2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups26.txt 

47 

2  H  7 

2.25  H7 

2  5  Hz 

3.1  Hz 

9.5  V 

1  V 

2.2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups27.txt 

48 

2  H  7 

2.25  Hz 

2  5  Hz 

3.1  Hz 

9.5  V 

1  V 

2.2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups28.txt 

49 

2  H  7 

2.25  Hz 

2.5  Hz 

3.1  Hz 

9.5  V 

1  V 

2.2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1500 

Down 

groups26.txt 

50 

2  H7 

2.25  Hz 

2.5  Hz 

3.1  Hz 

9.5  V 

1  V 

2.2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1400 

Down 

grou  ps26.txt 
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51 

2  Hz 

8.0  V 

5  cycles 

2.25  Hz 

1  V 

1  cycles 

2.5  Hz 

2.2  V 

1  cycles 

3,1  Hz 

8.0  V 

2  cycles 

1300 

Down 

groups29.txt 

52 

2  Hz 

2.25  Hz 

2.5  Hz 

3.1  Hz 

8.0  V 

1  V 

2.2  V 

8.0  V 

4  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups30.txt 

53 

2  Hz 

2.25  Hz 

2.5  Hz 

3.1  Hz 

8.0  V 

1  V 

2.2  V 

8.0  V 

3  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups31  txt 

54 

1.6  Hz 

2.25  Hz 

2.5  Hz 

3.1  Hz 

8.0  V 

1  V 

2.2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups32.txt 

55 

1.6  Hz 

2.25  Hz 

2.5  Hz 

3.1  Hz 

8.0  V 

1  V 

22  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups33.txt 

56 

2  Hz 

2.25  Hz 

2.5  Hz 

3.1  Hz 

9  5  V 

l  V 

2.2  V 

8.0  V 

5  cycles 

1  cycles 

l  cycles 

2  cycles 

1300 

Down 

groups26.txt 

57 

2  Hz 

2.25  Hz 

2.5  Hz 

3  1  Hz 

9.5  V 

1  V 

2.2  V 

8.0  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1 100 

Down 

groups26.txt 

58 

2  Hz 

2.25  Hz 

2  5  Hz 

3.1  Hz 

8  V 

1  V 

22  V 

6.5  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Dow  n 

groups34.txt 

59 

2  Hz 

2.25  Hz 

2.5  Hz 

3.1  Hz 

8  V 

0.5  V 

1.1  V 

6.5  V 

5  cycles 

1  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups35.txt 

60 

2  Hz 

2.25  Hz 

2.5  Hz 

3  1  Hz 

8  V 

I  V 

11  V 

6.5  V 

5  cycles 

2  cycles 

1  cycles 

2  cycles 

1300 

Down 

groups36.txt 

61 

2  l  lz 

2.25  Hz 

2.5  Hz 

3  1  Hz 

8  V 

1  V 

1.1  V 

6.5  V 

5  cycles 

2  cycles 

1  cycles 

2  cycles 

1500 

Down 

groups36.txt 

*  Run  numbers  omitted  were  either  runs  where  zeros  were  taken  or  bad  runs. 
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Table  12.  Phase  II  Wave  Run  Matrix 


Run 

Number** 

Wave  sequence 

file 

Blower 

RPM 

Lips 

Spectra  Embedded  Into 

GLRP 

Measurement 

05 

groups29.txt 

1300 

Down 

N/A,  repeat  of  run  5 1  from  Aug  07 

No 

06 

groups36.txt 

1500 

Down 

N/A,  repeat  of  run  61  from  Aug  07 

No 

07 

groups09.txt 

1300 

Down 

N/A,  repeat  of  run  1 5  from  Aug  07 

Yes 

08 

groups25.txt 

1300 

Down 

N/A,  repeat  of  run  43  from  Aug  07 

Yes 

09 

BR15130S.gen 

1440 

Down 

Bretsehneider  SS8,  >.=30, 

Hs=  15.0  in.,  Tm=  3.0s 

Yes 

10 

groups37.txt 

1300 

Down 

groups09.txt  at  end  of 

Bretsehneider  SS8,  >.=30 

Yes 

1  1 

groups38.txt 

1300 

Down 

groups09.txt  in  middle  of 

Bretsehneider  SS8,  >.=30 

Yes 

14 

groups39.txt 

1300 

Down 

groups25.txt  in  middle  of 

Bretsehneider  SS8,  >.=30 

Yes 

15 

ray01.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  o=  0.2,  3  V  max 

Yes 

16 

ray02.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  cr=  0.2,  6  V  max 

Yes 

17 

ray03.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  a=  0.2,  9  V  max 

Yes 

18 

ray04.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  a-  0.25,  3  V  max 

Yes 

19 

ray05.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  o=  0.25,  6  V  max 

Yes 

20 

ray06.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  o=  0.25,  9  V  max 

Yes 

21 

groups40.txt 

1300 

Down 

groups36.txt  in  middle  of 

Bretsehneider  SS8,  >.=30 

Yes 

22 

groups41.txt 

1300 

Down 

groups29.txt  in  middle  of 

Bretsehneider  SS8,  X  30 

Yes 

25*** 

Camlsr30.gen 

1600 

Down 

Hurrieane  Camille  speetra,  )i=30, 

Hs=  16.0  in.,  Tm=  2.45s 

Yes 

26 

Groups42.txt 

1600 

Down 

groups36.txt  in  middle  of 

Hurrieane  Camille  speetra,  X=30 

Yes 
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27 

groups42.txt 

1300 

Down 

groups36.txt  in  middle  of 

Hurricane  Camille  spectra,  X=30 

Yes 

28 

groups43.txt 

1300 

Down 

groups09.txt  in  middle  of 

Hurricane  Camille  spectra,  X-30 

Yes 

31*** 

rayl5.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  cr=  0.5,  9  V  max 

Yes 

32 

rayl3.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  cr=  0.5,  3  V  max 

Yes 

33 

rayl0.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  o=  0.4,  3  V  max 

Yes 

34 

ray07.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  cr  0.3,  3  V  max 

Yes 

35 

ray  1  l.txt 

1300 

Down 

N/A,  Rayleigh  distribution 

superposition,  o=  0.4,  6  V  max 

^  es 

36 

sincOI  .txt 

1300 

Down 

N/A,  Sine  wave,  frcq-0.3  Hz,  3  V 

max 

^  es 

37 

sinc04.txt 

1300 

Down 

N/A,  Sine  wave,  freq-0.5  Hz,  3  V 

max 

Yes 

38 

sinc07.txt 

1300 

Down 

N/A,  Sine  wave,  frcq  0.7  Hz,  3  V 

max 

Yes 

39 

sinc03.txt 

1300 

Dow  n 

N/A,  Sine  wave,  freq  0.3  Hz,  9  V 

max 

Yes 

40 

epszeaml.gen 

1600 

Down 

Hurricane  Camille  spectra,  X=46.6, 

Hs=  10.30  in  ,  Tm  1  96s 

Yes 

41 

groups44.txt 

1300 

Down 

groups36.txt  in  middle  of 

Hurricane  Camille  spectra,  X=46  6 

^  es 

42 

groups20.txt 

1300 

Up 

N/A,  repeat  of  run  35  from  Aug  07 

Yes 

43 

groups45.txt 

1 300 

Up 

groups20.txt  in  middle  of 

Hurricane  Camille  spectra,  X=46  6 

Yes 

**Run  numbers  omitted  were  either  runs  where  zeros  were  taken,  where  additional  repeat 
runs  were  performed,  or  were  bad  runs. 

***Runs  24  and  30  were  repeat  runs  of  14  and  28,  where  additional  GL.RP  measurements 
were  taken. 
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Appendix  B:  Senix  Results —  Regular  Waves 

Senix  probe  results  from  the  wave-packet  method  and  the  finite-wave  linear 
superposition  method  are  shown  in  this  Appendix.  Wave  group  characteristics  are  also  shown. 

Wave-Packet  Method 
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sec  sec 

Probe  6  Probe  3 
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Probe  4  Probe  1 


Figure  25.  Phase  II,  Run  36  Time  scries  data  from  sonic  probes  8,  6,  5,  4,  3,  and  1  with  the 

sinusoidal  wave-maker  signal  at  0.3  Hz 
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Figure  26.  Phase  II  Run  15  Time  series  data  from  sonic  probes  8,  6,  5,  4,  3,  and  1 

for  wave-packet  generation 
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Effect  of  Input  Amplitude 


Figure  27.  Input  signals  for  the  wave-maker.  Maximum  amplitudes  are  3,  6  and  9  Volts 
for  the  blue,  green  and  red  lines,  respectively. 
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Figure  28.  Phase  II,  Run  15  Time  series  data  from  sonic  probes  8,  6,  5,  4,  3,  and  1.  The 
probe  signals  were  multiplied  by  3, 1.5  and  1  for  the  blue,  green  and  red  lines,  respectively. 
Maximum  amplitudes  of  wave-maker  input  signals  are  3,  6  and  9  Volts  for  the  blue,  green 

and  red  lines,  respectively 
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wave  amp  (cm)  C  wave  amp  (cm) 


Probe  3 


re  29.  Phase  II,  Run  15  wave  amplitude  time  series  for  probe  3 


Probe  1 


Figure  30.  Phase  II,  Run  15  wave  amplitude  time  series  for  probe  1 
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Finite-Wave  Linear  Superposition  Method 


Probe  8 


Figure  31.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach).  Phase  I,  Run  15 
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Figure  32.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  I,  Run  35 


B-8 


Height  (cm)  Height  (cm) 


Probe  3 


Figure  32  (continued) 
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Figure  33.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  I,  Run  43 
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Figure  33  (continued) 
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Figure  34.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  I,  Run  51 
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Figure  35.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  I,  Run  61 
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Figure  36.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach).  Phase  II,  Run  5 
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Figure  37.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  II,  Run  6 
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Figure  37  (continued) 
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Figure  38.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach).  Phase  II,  Run  7 
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Figure  39.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach).  Phase  II,  Run  8 
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Figure  39  (continued) 
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Figure  40.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach).  Phase  II,  Run  42 
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Figure  41.  Wave  characteristics  at  sonic  probe  3,  Phase  I,  Run  15 
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Figure  42.  Wave  characteristics  at  sonic  probe  4,  Phase  I,  Run  15 
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Figure  43.  Wave  characteristics  at  sonic  probe  3,  Phase  I,  Run  35 
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Figure  44.  Wave  characteristics  at  sonic  probe  4,  Phase  I,  Run  35 
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Figure  45.  Wave  characteristics  at  sonic  probe  3,  Phase  I,  Run  43 
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Figure  46.  Wave  characteristics  at  sonic  probe  4,  Phase  I.  Run  43 
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Figure  47.  Wave  characteristics  at  sonic  probe  3,  Phase  I,  Run  51 
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Figure  48.  Wave  characteristics  at  sonic  probe  4,  Phase  I,  Run  51 
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Figure  49.  Wave  characteristics  at  sonic  probe  3,  Phase  I,  Run  61 
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Figure  50.  Wave  characteristics  at  sonic  probe  4,  Phase  1.  Run  61 
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Figure  51.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  5 
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Figure  52.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  5 
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Figure  53.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  6 
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Figure  54.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  6 
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Figure  55.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  7 
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Figure  56.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  7 
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Figure  57.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  8 
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Figure  58.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  8 
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Figure  59.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  42 
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Figure  60.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  42 
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Appendix  C:  Senix  Results —  Irregular  Waves 

Results  from  the  Phase  II  testing,  with  embedded  wave  groups  in  irregular  waves,  are 
shown  in  this  Appendix.  The  results  include  Senix  probe  time-histories  and  analysis  of  the  wave 
group  characteristics. 
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Figure  61.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  II,  Run  1 1 
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Figure  61  (continued) 
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Figure  62.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  II,  Run  14 
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Figure  63.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  II,  Run  21 
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Figure  64.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  11,  Run  22 
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Figure  64  (continued) 
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Figure  65.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach).  Phase  II,  Run  24 
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Figure  66.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to 

(closest  to  the  beach),  Phase  II,  Run  26 
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Figure  66  (continued) 
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Figure  67.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  II,  Run  27 
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Figure  67  (continued) 
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Figure  68.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  II,  Run  28 
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Figure  68  (continued) 
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Figure  69.  Measured  nave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  II,  Run  30 
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Figure  70.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach).  Phase  II,  Run  41 
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Figure  70  (continued) 
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Figure  71.  Measured  wave  time-history  of  sonic  probes  8  (closest  to  wave-maker)  to  1 

(closest  to  the  beach),  Phase  II,  Run  43 
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Figure  71  (continued) 
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Figure  72.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  1 1 
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Figure  73.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  1 1 
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Figure  74.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  14 
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Figure  75.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  14 
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Figure  76.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  21 
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Figure  77.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  21 
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Figure  78.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  22 
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Figure  79.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  22 
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Figure  80.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  24 
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Figure  81.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  24 
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Figure  82.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  26 
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Figure  83.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  26 
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Figure  84.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  27 
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Figure  85.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  27 
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Figure  86.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  28 
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Run  28  Probe  4  - Wave  1  - Wave  2 - Wave  3 
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Figure  87.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  28 
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Figure  88.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  30 
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Run  30  Probe  4  - Wave  1  - Wave  2 - Wave  3 
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Figure  89.  Wave  characteristics  at  sonic  probe  4,  Phase  11,  Run  30 
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Run  41  Probe  3 - Wave  1  - Waw  2 - Ways  3 
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Figure  90.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  41 
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Figure  91.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  41 
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- Run  43  Probe  3 - Wave  1  - Wave  2 - Wave  3 
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Figure  92.  Wave  characteristics  at  sonic  probe  3,  Phase  II,  Run  43 
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- Run  43  Probe  4 - Wave  1  - Wave  2 - Wave  3 
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Figure  93.  Wave  characteristics  at  sonic  probe  4,  Phase  II,  Run  43 
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Figure  94.  Spectral  Analysis  of  Phase  II,  Run  1 1-  Bretschneider  SS8, 1=30,  lls=  38.1  cm 
(15.0  in.),  Tm=  3.0s  with  large  grouped  waves  embedded 
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Figure  95.  Spectral  Analysis  of  Phase  II,  Run  26-  Hurricane  Camille,  X=30,  Hs=  40.64  cm 
(16.0  in.),  Tm=  2.45s  with  wave  group  embedded 


C-47 


Figure  96.  Spectral  Analysis  of  Phase  II,  Run  40-  Hurricane  Camille,  X=46.6,  Hs=  26.16  cm 
(10.3  in.),  Tm=  1.96s  with  wave  group  embedded 
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Appendix  D:  GLRP  Results 

Results  from  GLRP  measurements  and  comparisons  to  Senix  probe  data  are  shown  in  this 
Appendix.  Three-dimensional  representations  of  GLRP  data  are  presented  every  third  figure,  in 
Figures  97-121.  These  representations  show  several  snapshots  of  the  wave  group  evolution  in 
the  measurement  region.  The  first  frame  shows  the  water  surface  before  the  wave  group  enters 
the  measurement  region.  The  second  frame  shows  the  beginning  of  the  wave  group.  In  the  third 
and  fourth  frames,  the  maximum  and  minimum  wave  heights  encountered  in  the  measurement 
are  shown.  The  final  frame  shows  the  ealm  water  surface  at  the  end  of  the  run. 
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Figure  97.  Three-dimensional  wave  surface  plots  for  Phase  II,  Run  7 
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GLRP  and  Sonic  Probe  5  data  for  Diode  at  position  X  =  3462.37mm,  Y  =  1304  03mm 


Figure  99.  GLRP  single  point  and  Sonic  data  comparison  for  Phase  II,  Run  7 


D-3 


(0,5)7 


QCRP  R  m  f<v  fl .  A I  CMimi 
Mn«km  W  rr«  E 

(■  II  0«7  •  « 


04.  R^  Rm»(i  f**  Rur  t  AICir**n« 
Mri  *%  M  iw 
I  ■  W  70*i*< 


y.  '0 


to  '*  :o  :*g  so  » 


#  :»  :c  'i  i  o 


k  r. 


tlW  NMuh^fiunl.M  Cmmhi 
A**#  W»v*  Qr»wpP»ti*« 

(•80  CMik 


'*>  .v  .*  t»  .‘C  *  0  1  k  •!*  :o  .*  w  .« 


Figure  100.  Three-dimensional  wave  surface  plots  for  Phase  11,  Run  8 
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'igure  101.  GLRP  Single  Point  Time  History  for  Phase  II.  Run  8 

GLRP  and  Sonic  Probe  5  data  for  Diode  at  position  X  =  3465.05mm.  Y  =  1303  28mm 


Figure  102.  GLRP  single  point  and  Sonic  data  comparison  for  Phase  II,  Run  8 
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Figure  103.  Three-dimensional  wave  surface  plots  for  Phase  II,  Run  11 
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Figure  104.  GLRP  Single  Point  Time  History  for  Phase  II,  Run  1 1 


GLRP  and  Sonic  Probe  5  data  for  Diode  at  position  X  =  3460  85mm,  Y  =  1386  44mm 


Figure  105.  GLRP  single  point  and  Sonic  data  comparison  for  Phase  II,  Run  1 1 
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Figure  106.  Three-dimensional  wave  surface  plots  for  Phase  II,  Run  21 
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GLRP  and  Sonic  Probe  5  data  for  Diode  at  position  X  =  3548  25mm,  Y  =  1385  74mm 


Figure  108.  GLRP  single  point  and  Sonic  data  comparison  for  Phase  II,  Run  21 
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Figure  109.  Three-dimensional  wave  surface  plots  for  Phase  11,  Run  22 
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re  110.  GLRP  Single  Point  Time  History  for  Phase  II.  Run  22 


GLRP  and  Sonic  Probe  5  data  for  Diode  at  position  X  =  3462  63mm,  Y  =  1386  04mm 


Figure  111.  GLRP  single  point  and  Sonic  data  comparison  for  Phase  II,  Run  22 
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Figure  112.  Three-dimensional  wave  surface  plots  for  Phase  II,  Run  24 
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Figure  1 13.  GLRP  Single  Point  Time  History  for  Phase  11,  Run  24 


GLRP  and  Sonic  Probe  5  data  for  Diode  at  position  X  =  3463.06mm,  Y  =  1386  74mm 


Figure  1 14.  GLRP  single  point  and  Sonic  data  comparison  for  Phase  II,  Run  24 
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Figure  1 15.  Three-dimensional  wave  surface  plots  for  Phase  II,  Run  26 
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GLRP  and  Sonic  Probe  5  data  for  Diode  at  position  X  =  3462.00mm,  Y  =  1303  63mm 
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Figure  118.  Three-dimensional  wave  surface  plots  for  Phase  II,  Run  27 
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Z  at  Diode  X  =  3463.17mm,  Y  =  1386  49mm 


Figure  119.  GLRP  Single  Point  Time  History  for  Phase  II,  Run  27 

GLRP  and  Sonic  Probe  5  data  for  Diode  at  position  X  =  3463  17mm,  Y  =  1386  49mm 


1  \ 

1 - 1 - » - 

1  1 

L— -  *  ~ 

GLRP  point 
Sonic  5  data 

Figure  120.  GLRP  single  point  and  Sonic  data  comparison  for  Phase  II,  Run  27 
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Figure  121.  Three-dimensional  wave  surface  plots  for  Phase  II,  Run  30 
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Z  at  Diode  X  =  3462  98mm,  Y  =  1385  84mm 


GLRP  and  Sonic  Probe  5  data  for  Diode  at  position  X  =  3462.98mm,  Y  =  1385  84mm 


Figure  123.  GLRP  single  point  and  Sonic  data  comparison  for  Phase  II,  Run  30 
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